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REACTOR SAFETY 


Highlights of the Aix Conference 


on Fast Reactor Safety 


By C. E. Dickerman 


A general spirit that safety problems offer no basic 
obstacle to development of sodium-cooled fast reac- 
tors permeated the recent International Conference 
on Fast Reactor Safety held at Aix-en-Provence, 
France.! Speaker after speaker called for vigorous 
moves to develop large fast power reactors. 
W. McCarthy, of the Power Reactor Development 
Company, Fermi reactor, admonished that the fast 
breeder should become the standard reactor of the 
future: “the reactor our children will know.” 

Unlike the 1965 fast reactor safety conference at 
Argonne,” when reactivity coefficients were given the 
most technical attention, this September 1967 Aix 
conference was noted by the vigorous attention to 
the models for describing and calculating transient 
sodium boiling and expulsion by vaporization. The 
interest culminated in an informal meeting of the 
sodium investigators attending, who recommended 
that a conference on sodium boiling should be held 
as a follow-up. These developments will be reported 
in detail in a subsequent article. 

Problems of fuel meltdown, sodium-vapor explo- 
sion, and malfunction detection are recognized but 
now are considered comparable to those of high- 
power-density thermal reactors. D. Okrent of Argonne 
cautioned against overlooking safety problems and 
observed that experience has shown that the unex- 
pected often does occur. G. Denielou of the French 
Commissariat 4 l’Energie Atomique summed up his 
feelings by noting that problems discussed seem to 
be essentially the same as for thermal reactors. 
F. Farmer of the UKAEA continued in a parallel 
vein by expressing the U. K. attitude that the fast 
reactor is proving to be “just another highly rated 
reactor” system. 

Topics discussed include the physics of both local 
and generalized accidents, neutronics, critical ex- 
periments, sodium boiling, loss of coolant, reactor 
operation, and accident detection. The low ratio of 
steam papers (4) to sodium papers (63) reflected 
the actual emphasis of the discussions. Of the na- 


tional programs represented at the conference, only 
the German program gave comparable emphasis to 
both steam-cooled reactors and sodium-cooled re- 
actors. 


Steam-Cooled Reactors, Briefly 


The steam-cooled system has evolved into a direct 
cycle, with the reactor supplying superheated steam 
to the turbine; thus it more closely resembles gas- 
cooled-reactor designs than sodium-cooled fast re- 
actors. 

A paper from BelgoNucléaire focused on safety 
aspects of direct steam cycles.° Two German papers!»® 
and a U. S. paper® presented detailed considerations 
of equipment and control characteristics of the 
coupled reactor—turbine system. 

The conference papers showed clearly the relative 
stage of development of steam-cooled fast reactors— 
systems now in the design-study stages— and liquid- 
metal-cooled fast reactors, which first moved to the 
operating stage a generation ago with Experimental 
Breeder Reactor I (EBR-I). No significant safety 
advantages of the steam systems were apparent. 


Liquid-Metal Systems 


It is apparent from the desire for details, both in 
experimental studies and in theoretical analyses, that 
the technology of liquid-metal-cooled fast breeder 
reactors is much more advanced than is the technol- 
ogy of steam-cooled breeders. In the liquid-metal 
systems, the central questions relate to meltdown, 
core compaction, and sodium behavior as they are 
characterized by experiment, analysis, and operating 
fast reactor experience. 


MELTDOWN EXPERIMENTS 


Three papers on experimental studies of core 
behavior under accident conditions were presented. 
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First was an Argonne report on its Transient Reactor 
Test (TREAT) reactor-transient studies, which fea- 
tured descriptions and analyses of recent meltdown 
experiments in sodium loops.* The loop tests were 
performed on metallic, EBR-II driver fuel pins, run 
both as single pins and in small seven-pin clusters. 
Fuel-failure thresholds were reviewed in terms of a 
basic single-pin model. Rapid sodium expulsion, sharp 
pressure spikes, and extensive secondary movement 
of fuel and coolant were reported. TREAT was de- 
signed specifically for short transients, with durations 
of the order of a few seconds. 

The second paper, from Cadarache, described the 
CEA program under way to set up a sodium loop for 
slow (up to 10 min) oxide-fuel meltdown experiments 
on the Cabri swimming-pool reactor.® The test loop 
being developed is designed for experiments on single 
fuel pins, with a capability for expansion to three-pin 
clusters by means of minor changes. The goal is 
meltdown experiments on full subassemblies, and a 
reference case of a Seven-subassembly cluster was 
given as the eventual loop design basis. 

Next was a BelgoNucléaire report? that analyzed 
results of a power excursion on oxide fuelin a sodium 
loop. The analysis was performed as part of the 
BelgoNucleaire studies on their proposed transient 
test reactor Briseis. 


CORE-COMPACTION ACCIDENTS 


It has become traditional for fast reactor safety 
sessions to include papers on the calculation of 
nuclear explosions that result from superprompt- 
critical reactivity insertions brought about by melt- 
down and core compaction during a less severe 
accident. A German paper!’ presented results from 
calculations on a mock-up of the core pin geometry 
for the core disassembly. Their geometry model 
was chosen as an improvement over the usual hy- 
pothetical fluid assumed for such computations. The 
paper by Lord'! considered effects of geometry 
changes, but with a conventional core material de- 
scription. In general, the calculations indicate that 
a value of about 500 lb “TNT equivalent” is a rea- 
sonable basis for containment design of prototype- 
size fast reactors. 


What Core Behavior? Core behavior leading to the 
explosion phase was not detailed in these explosion- 
calculation papers. If a nuclear explosion is assumed 
to result from an initiating, less-severe accident in- 
volving all (or a large part) of a core, it is necessary 
to examine the course of the initiating accident, in- 
cluding the degree of coherence in fuel failure, of 
fuel and coolant movement, and, finally, of possible 
secondary collapse of fuel remains; if the accident 
is assumed to result from a severe localized acci- 
dent (such as subassembly blockage and meltdown 
at power), it is necessary to consider the severity 
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and speed of propagation of damage and degree of | 


coherence of secondary movements. Thus it is clear — 
to this reviewer, at least—that the biggest current 
problem in the maximum-accident nuclear explo- 
sion calculations is the description of core behavior 
leading to the maximum accident. A related problem 
is the adequacy of the physical description of the 
core during the explosion. 


SODIUM VAPORIZATION 


Special attention was given to the problem of rapid 
energy transfer from hot finely divided fuel fragments 
to sodium. The resulting work done by the expansion 
of sodium vapor is analogous to a steam explosion 
such as those generated in the BORAX and SPERT 
experimental water-reactor destructive tests. Farmer 
(United Kingdom) observed in his panel talk (at the 
conclusion of the conference) that the thermal energy 
remaining in a fast reactor core after a maximum- 
accident explosion is very large and that it is possi- 
ble to calculate a sodium-vapor explosion as severe 
as the (assumed) previous nuclear excursion. Hence 
the United Kingdom considered it as important to 
design containment for a sodium-vapor explosion as 
to design for a core vaporization explosion. 


Models, MacFarlane'* described calculations with 
a coupled code system that combines transient heat 
transfer, two-phase coolant hydrodynamics using an 
annular flow model, and neutron kinetics. Although the 
paper was presented to display a significant advance 
in the development of accident analysis techniques, 
the conference interest in details of transient sodium 
boiling was so strong that the vigorous discussion 
following the paper was almost exclusively on the 
two-phase hydrodynamics module of the coupled code 
system! 

A contrasting model, of sodium slug expulsion, was 
incorporated into a Karlsruhe sodium-voiding code, 
the BLOW code, which was described in the paper by 
Peppler.'? The BLOW model was developed on the 
basis of analyses of experiments performed by Grass 
and coworkers at Ispra!4 and was also related to re- 
sults obtained by LeGonidec under CEA auspices.'* 

The effects of two-phase compressibility phenomena 
on transient sodium behavior and analyses of data 
were discussed by Fauske.'® Liquid-metal superheat 
values reported! ranged up to 800°C, with widely 
varying but lower values reported by French,'° 
Argonne,'’ Atomics International,'® and Casaccia!® 
workers. However, Pinchera’® pointed out that under 
realistic core conditions the effective value of super- 
heat might be rather low. 

The theory of transient vaporization and condensa- 
tion of sodium was discussed in some detail by 
Peppler,'® Judd,”’? and Kosky.”! For the models he 
used, Judd?’ found that dynamics of void growth were 
relatively insensitive to changes in geometry ranging 
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from a single spherical bubble to many spherical 
bubbles. 

The general problem of sodium behavior during a 
core accident is still not well defined. Because of 
the scope and complexity of the subject, the Aix 
sodium papers will be discussed in greater detail 
in separate reviews on reactor safety implications 
of superheat” and on coolant-ejection mechanisms. ”? 


Sodium Release. Consequences of sodium release 
during an accident also were considered. A paper 
by Morewitz™ reported data and analyses of aerosol 
generation and behavior during a sodium fire. Their 
results showed that the presence of sodium oxide 
particles produced a significant scavenging action 
on the fuel aerosols. Deschamps and Rygaert” de- 
scribed related experiments on iodine distribution 
following a simulated sodium release. Containment 
of plutonium aerosols was discussed by Heller. 7° 


OPERATING EX PERIENCE 


A high point of the conference was the session 
on actual operating experience with the Dounreay, 
Rapsodie, BR-5, and Fermi fast reactors. 


Dounreay. Phillips*’ reviewed the history of the 
Dounreay reactor, which was shut down by aNaK leak 
shortly before the conference. The paper was essen- 
tially a discussion of the well-known use ofthis facil- 
ity as a fast-flux testing reactor. There was little 
more than an acknowledgment concerning the NaK 
leak; however, the leak has recently been located in 
one of 24 links between the primary coolant circuit 
and the reactor vessel.”® Repair of the leak is ex- 
pected to require ‘some months,” apparently as a 
result of accessibility and radiation problems. 


Rapsodic. Gajac and Reynes”® told of problems 
encountered in the commissioning and startup of the 
French Rapsodie reactor, including thermal-stress 
difficulties in the original design and an inadvertent 
release of sodium. The actual reactor operation was 
reported to be smooth; in fact, current French opti- 
mism on the state of the fast-breeder-reactor art 
was such that the prototype Phenix reactor has been 
designed without a sealed containment building. *” 


BR-5. Orlov*® told of the variety of difficulties 
met and overcome in the operation of the Russian 
BR-5 reactor. The report presented new informa- 
tion on BR-5 experience, such as extensive oxide 
plugging of the core following exposure of a drained 
core to air, gas entrapment in the primary system, 
and leaks in the liquid-metal heat exchanger. 


Fermi. Most interest in the operating experience 
was generated by the review of the partial meltdown 
in the Fermi reactor. McCarthy” described the 
sequence of events in the incident, which occurred 


a year before the meeting, and announced that the 
cause— an unidentified object in the inlet plenum— 
had been found just prior to the conference. As 
follow-up papers to the Fermi review, there were 
analyses on thermocouple response in Fermi to core 
abnormalities®® and analyses of instrumentation that 
could have detected the Fermi meltdown (which oc- 
curred on startup) in its early stages. *4 The Fermi 
operators have been quite timely and candid in their 
reports on the incident. A somewhat condensed ver- 
sion of their Aix report®” has been published. *® 
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DESIGN AND CONSTRUCTION PRACTICE 


Further Developments in Containment-Building 
Electrical Penetrations 
By Fred Verber 


New information about the electrical penetration sys- 
tems in the Fermi,! Peach Bottom,” and Dresden 2 
(Ref. 3) nuclear power stations supplements that pre- 
sented previously in this journal.‘ Specifically, this 
article considers some electrical penetration prob- 
lems and solutions experienced during plant opera- 
tions at Fermi and at Peach Bottom, describes fire- 
protection modifications to the penetrations and tothe 
electrical design at Peach Bottom following the oc- 
currence of a fire, and presents recent developments 
in the design of the Dresden 2 electrical penetration 
system and assemblies. 


Fermi 


Figures 1 to 3 are views of Fermi electrical pene- 
trations from inside and outside the reactor building. 
This containment vessel is designed for an internal 





Fig. 1 Some Fermi electrical penetrations and spares are 
viewed from outside reactor containment building. 








? 


Fig. 2 Packing nut being installed (right) on outside end of 
Fermi penetration nozzle.’ 


| 





Fig. 3 Inner end of two Fermi penetrations, showing both 
tubing connections for continuously monitoring for leakage 
and extrusion of backup sealant! 
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pressure of 32 psig. The cluster of penetrations shown 
in Fig. 1 is part of a total of 517 penetrations (ap- 
proximately 325 in use) grouped above the operating 
floor in a sector of the building that is 35 ft long and 
15 ft high. The penetration nozzles are spaced ap- 
proximately 10 in. on centers. These nozzles carry 
cables for 4800-, 480-, and 120-volt and lower voltage 
circuits; cable diameters range from '/, in. to 2 in. 


Figure 2 shows two stages of the installation of a 
packing nut that seals the cable at the outside end of a 
penetration nozzle. Figure 3 shows an interior view of 
two penetrations. Note the extrusion of the backup 
sealant (trade name, Weatherban) on the penetration; 
also note the tubing connections for continuously mon- 
itoring the penetrations for leakage. 


Cold-Weather Effects. As previously reported,‘ 
the application of Weatherban, a butyl- rubber resilient 
sealant, to the neoprene packaging rings and other 
parts of the penetration assembly contributed very 
significantly to the integrity of the penetration seals. 
But during extremely cold weather, the outer seal 
apparently loses some of its resiliency, and occa- 
sional leaks have resulted.' These seals were taken 
apart and reassembled using new packing and caulking 
with Weatherban. The header system, which supplies 
5 psi air pressure continuously to each penetration, 
is equipped with a flowmeter that is checked daily for 
indications of increased leakage. Each month the 
penetrations are tested at 32 psig. 


A leak-rate test of the electrical penetrations was 
conducted at 32 psig on Nov. 6, 1967, and the total 
leakage was 551 scf/day. After repairs were made on 
the penetrations with excessive leakage, a leak-rate 
test was performed later in the month, and the total 
leakage was 291 scf/day. This was followed by ad- 
ditional repairs. No repairs had been required during 
the preceding 7 months. During the year from Feb- 
ruary 1967 to February 1968, repairs were made on 
11 of the 325 penetrations in use.» 


Peach Bottom 


Containment of the reactor and auxiliaries of the 
Peach Bottom Atomic Power Plant is provided by a 
steel vessel designed to withstand an internal pres- 
sure of 8 psig. 

During plant construction a fire occurred in one of 
the U-tube cable penetrations inside the reactor con- 
tainment building.” It started when hot weld metal 
dropped onto the cable penetrations as a platform was 
being constructed above them. Because of the fire, 
several modifications have been made to improve fire 
protection and plant safety. 

About 1100 cables were damaged in the fire. Of 
these, the multiconductor control and instrument 
cables were insulated with polyethylene and jacketed 
with polyvinyl chloride, the power cables were insu- 
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lated with butyl rubber and jacketed with polyvinyl - 


chloride, and the thermocouple extension leads were 
jacketed with polyethylene. Tests showed that poly- 
vinyl chloride will burn but will not support combus- 
tion at moderate temperatures. However, the poly- 
ethylene and the permagum, silicone, or thiokol filler 
materials in the cables are combustible. 
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Fig. 4 Peach Bottom thermocouple-cable gas stop shows 
(W~ vermiculite fire-proofing malerial added at ends of U- 
lube pene tralions and (b) « boxy resin that replaced grout Jor 
sealing individual wires of unfilled mulliconductor cables. 


FIRE-PROTECTION IMPROVEMENTS 

Modifications made to improve plant safety include 
the following: 

e For reliability, certain power and control cables 
were rerouted to segregate the circuits of essential 
redundant equipment. 

¢ Horizontal runs of cable trays were covered for 
protection against falling objects and to prevent the 
spread of fire from one tray to another. 

¢ Metal covers and fire stops were installed in long 
vertical runs of cable trays. 

e Vermiculite, a fireproofing material, was added 
to both ends of the U-tube penetrations as shown in 
Fig. 4. 

e Smoke detectors were installed in areas where 
electrical cables and equipment are concentrated, in- 
cluding the 480- and 2400-volt switchgear rooms and 
motor— generator sets. 

e The CO, fire-extinguishing system was extended 
within the containment building to include areas where 
electrical cables are heavily concentrated. 
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Several other changes were made to improve leak- 
tightness of the cable penetrations: 

e Epoxy resin was substituted for grout for sealing 
individual wires of unfilled multiconductor thermo- 
couple cables after the outer jacket was removed 
(Fig. 4). 

e Solid connectors were installed on single- 
conductor 600-volt cables to prevent leakage between 
strands of the conductor; the connector was covered 
by a heat-shrinking Raychem sealing sleeve (Fig. 5). 
This connector system was found to be very effective. 

e Where leakage through cables was detected during 
tests, gas stops, as shown in Fig. 6, were added. 

Individual leak-rate tests of all new penetrations 
showed the leakage to be insignificant. An integrated 
containmen*-building leak-rate test at 8 psi showed 
that the leak rate was less thanthe specified 0.2%/day 
and that the containment integrity of the building had 
not been impaired by the fire. 

These modifications represent significant design 
improvements and safeguards, particularly in respect 
to fire protection of electrical equipment and wiring. 
Many of these features are applicable to— and should 
significantly improve equipment and wiring installa- 
tions of— power-generating stations and industrial 
plants as well as reactor containment buildings. 
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Fig. 5 Solid connectors installed on single-conductor un- 
filled cables prevent leakage between strands of conductor 
at Peach Bottom. Raychem heat-shrinking sealing sleeves 
cover connectors. 
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al Peach Bottom where leakage had been detected through 
cable.” 


Dresden 2 


The design of the electrical penetration system for 
the Dresden 2 primary containment system has ad- 
vanced since the earlier description. The primary 
containment design employs a pressure-suppression 
system that includes a steel pressure-vessel dry well 
containing the reactor vessel and recirculation sys- 
tem. The design internal pressure of the dry well is 
62 psig. A reactor building provides secondary con- 
tainment; it encloses the primary containment system 
and reactor servicing and auxiliary equipment. The 
design internal pressure of the reactor building is 
0.25 psig. The electrical cables entering the building 
will pass through ducts sealed in the wall; soft sealing 
materials will be used around these cables. 

As shown in Table 1, there are a total of 8376 con- 
ductors penetrating the dry-well structure in the 
present development of the Dresden 2 design; of these, 
5568 are for control-rod-position indication, 448 for 
reactor gamma and neutron monitoring, 2348 for 120- 
and 440-volt power, control, and instrumentation cir- 
cuits, and 12 for high-voltage (4160-volt) power cir- 
cuits. Approximately 18% of the conductors are spares. 

At an early stage of the design, a total of 63 elec- 
trical penetrations were contemplated. As the table 
indicates, the number of penetrations has been re- 
duced to 20. In addition, seven spare penetration 
sleeves are provided (Figs. 7 and 8) in the dry-well 
structure. 
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Fig. 7 Cutaway of Dresden 2 
other service penetrations are also shown. 
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Table 1 DRESDEN 2 ELECTRICAL PENETRATION SYSTEM? 

















Total No. of 
Type of penetration number of spare 
assembly Quantity Service Range of conductor sizes conductors conductors 
Low-voltage power 6 120-, 240-, and 480-volt #14 to #0 wires, #20 2348 1006 
and control cable power, valves, limit thermocouple wire, 
penetration ewitches, thermocouples, 400-MCM bare 
cortrols ground wire 
Low-voltage power 6 Control- rod-position Thirty-two 27-conductor 5568 135 
and control cable indication #18 wires plus 1 pair of 
penetration Chromel— Alumel thermo- 
couples per penetration 
High-voltage power t 4160-volt recirculating- Three 1000-MCM, single- 12 0 
cable penetration pump drive motors conductor cables per 
penetration 
Shielded-cable 4 Reactor neutron and Fifty-one type RG-108 A/U 448 96 
penetration gamma monitoring shielded pairs; 5 coaxial 
triple-shielded cables 
Totals 20 8376 1537 


Miscellaneous low - voltage 
power and control leads 













Control-rod-position indication 
Spare sleeve 


4160-volt recirculating- 
pump motor leads 


Control-rod-position 


indication —-Q 


Spare sleeve 


Miscellaneous low-voltage 
power and control leads 


Welding receptacles and 
dry-well cooler blowers 


Spare sleeve 


Dry - well steel shell 


| ~_— 32.0 ft radius & 
q 38.25 ft radius 
D 
q 4 
D> v 
D Z > 
e: ™. 
> % » 
hy, Ql <s Or 
CP ‘ 4160-volt recirculatin 
Control-rod-position rare 4 pond motor leads ' 
indication en ae 


Thermocouples 


Neutron-detector circuits 
Spare sleeve 





Source and intermediate-range neutron- 
* Z d r-dr t 
Neutron-detector circuits L) etector-drive circuits 


Control-rod-position indication 


Fig. 8 Horizontal section through Dresden 2 dry-well structure (at level shown in Fig. 7) shows relative 
location of all 20 electrical penetration assemblies and seven spare sleeves.$ 


REACTOR AND FUEL-PROCESSING TECHNOLOGY, Vol. 11, No. 1, Winter 1967— 


1968 





10 DESIGN AND CONSTRUCTION PRACTICE 


TYPES OF PENETRATIONS 

A family of three commercially available standard 
types of electrical penetration assemblies, as devel- 
oped by General Electric Company for application in 
the nuclear power industry, will be used in the pri- 
mary containment structure for Dresden 2. The basic 
configuration of the three types of assemblies is es- 
sentially the same. Each penetration assembly is in- 
serted into a penetration nozzle (Fig. 9) thatis part of 
the containment dry well; double seals are formed by 
the two header plates of the assembly, thuscreating a 
pressure chamber that facilitates leak-rate testing 


3-in. - thick end shield Field weld before 




















Containment 


and continuous monitoring. Each penetration assem- 
bly will be equipped with a drain valve and pressure 
gauge, as indicated in Fig. 9. The pressure chambers 
will be purged with dry nitrogen gas at 15 psig for 
monitoring during normal plant operation. 

Each of the three types of penetration assemblies 
is designed to satisfy the needs of a particular group 
of electrical cables. The three types are categorized 
as follows: 


Low-voltage power and control cable penetration. 
This penetration (Fig. 10) is used for a-c circuits up 
to 600 volts and wire sizes from AWG #18 to 4/0. 
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Fig. 9 Typical Dresden 2 electrical-penetration-assembly canister is inserted into nozzle on containment 
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Low-voltage power and control cable penetration® at Dresden 2. 
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Fig. 12 Typical shielded-cable penetration at Dresden 2 handles reactor neutron-monitoring circuits.’ 


High-voltage power cable penetrations. This pene- 
tration (Fig. 11) accommodates circuit voltages up to 
5000 volts and cable sizes up to 1000 MCM. Non- 
magnetic stainless-steel header plates are used to 
eliminate the possibility of eddy-current heating 
(carbon-steel header plates are normally used in the 
other types of penetration assemblies). Three or six 
conductors per assembly is a typical conductor den- 
sity for this penetration. 


Shielded-cable penetration. Shielded signal cables 
are provided in this penetration (Fig. 12) to connect 
low-noise circuits between the reactor and control 
room. Coaxial connectors are mounted directly onthe 
header plates and are isolated from ground. Single-, 


double-, or triple-shielded coaxial cables are used 
with these connectors and are normally used for cir- 
cuits with a pulse output signal. Leads for reactor 
monitoring circuits with a d-c current signal consist 
of a twisted pair covered by a braid shield. 


Figure 7 shows the design arrangement of a group 
of electrical penetration assemblies and other service 
openings on the dry-well structure. Figure 8 shows 
all 27 electrical penetrations. 


Physical and operational characteristics of these 
three basic types of penetration assemblies are: 


e The penetration assembly has a maximum leak 
rate of 1.16 x 107° em*/sec when pressurized to 63 


REACTOR AND FUEL-PROCESSING TECHNOLOGY, Voi. 11, No. 1, Winter 1967-1968 








12 DESIGN AND CONSTRUCTION PRACTICE 


psig with dry helium in an ambient temperature of 
175° F. 

¢ The assembly will continuously withstand an in- 
ternal pressure of 125 psig during normal operation. 

¢ The assembly maintains containment integrity 
continuously with the following conditions inside the 
containment (dry well): 150° F, —2 to 2 psig, 20 to 
100% relative humidity, and 10 r/hr radiation dose 
(without shielding). 

¢ The assembly will maintain containment integrity 
for at least 2 hr under the following maximum condi- 
tions inside the containment dry well: 320° F, 125 
psig, and 100% relative humidity. 

¢ Each conductor of the assembly is cut and recon- 
nected in a solid mechanical splice at each header 
plate (best illustrated in Fig. 11). This eliminates 
leakage through the strands of the conductor. Splices 
and spaces surrounding the cables penetrating the 
header plate at each end of the assembly are buried in 
a potting compound, namely, a version of 3M Scotch- 
cast epoxy resin No. 8. The resin is applied under a 
vacuum to eliminate voids that otherwise would be 
formed by trapped bubbles. 

e The materials comprising the assembly are fire 
resistant and self-extinguishing in accordance with 
ASTM-D635. 

e The penetration assembly is designed and fabri- 
cated in accordance with ASVE Boiler and Pressure 
Code, Section III, Class B Vessel. 

e As indicated in Fig. 9, cable boxes are installed 
on both ends of the assembly to accommodate the ex- 
tension of cable leads in splices or in environmental 
resistant-type connectors. Figure 13 shows the fac- 
tory termination of penetration cables in 10 con- 
nectors installed in a cable box mounted on one end 
of a penetration assembly. 

¢ Protection against line-of-sight radiation is pro- 
vided by end shields consisting of 3-in.-thick carbon- 
steel disks mounted on the shield support structure 
(Fig. 9) installed on the inside ofthe dry well; in addi- 
tion to the disks, 2-in.-thick carbon-steel shielding is 
provided on the collars of the penetration assemblies. 


I thank E. L. Alexanderson (of Power Reactor De- 
velopment Company), R. F. Griffin (of Bechtel Corp,), 
and D. McDaniel and R. B. Hubbard (of General Elec- 
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Fig. 13 Factory assembly of one type of penetration as- 
sembly to be used in Dresden 2, Note termination of cables 
in 10 connectors in cable box. 
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Realistic Data Are Needed for Mechanical 
Analysis of Reactor Fuel Elements 


By Richard A. Valentin 


As reactor design has progressed into the develop- 
ment of high-power-density systems with greater 
burnup and thermal efficiency, problems in the me- 
chanical design of fuel elements have become more 
complex. One indication of this is the growing number 
of papers that show (1) application of rational methods 
of analysis to the theoretical prediction of thermal 
displacements, stresses, and strains, and (2) attempts 
to include radiation-induced creep and swelling effects 
in a more logical and less empirical manner. The 
analytical methods used in these studies have become 
more sophisticated and often include relatively recent 
developments from the mathematical theory of plas- 
ticity, viscoelasticity, and viscoplasticity. 

Rather than attempting to survey comprehensively 
the extensive literature on the mechanical analysis of 
reactor fuel elements, this article discusses a few 
recent papers to illustrate the salient features and 
typical difficulties of such analyses. Noteworthy are 
the many assumptions needed to produce tractable 
problem formulations. A sequence of closely interre- 
lated papers by Ma and his associates'~* is repre- 
sentative of this type of mechanical analysis, Of this 
group, we will consider “Mechanical Analysis for 
Solid Cylindrical Fuel Elements”! and “The Radiation 
Growth and Swelling Analysis for Cylindrical Fuel 
Elements During a Positive Reactor Period’ as 
typical examples. 


Time Ranges. A point that deserves particular 
emphasis is the time range over which the analysis 
is to be applicable. Perhaps no other factor has such 
a profound influence on the mathematical models that 
can be employed and the analytical methods available 
for problem solutions, Although no truly satisfactory 
classification of time ranges can be stated, four 
rough classes will be used for illustration: 

A. Very-short-term dynamic analysis of transient 
effects wherein one does not neglect inertia effects. 
In such problems, elastic wave propagation plays an 


important role in the behavior of the reactor. A 
typical example is the effect of wave propagation on 
the self-limiting behavior of pulsed fast reactors of 
the Godiva type (see Refs. 6 and 7). Time scales of 
this order, of course, preclude any interest in creep 
as the term is normally applied; however, it is 
conceivable that the safety analysis of such systems 
could include problems in the dynamic theory of 
plasticity or viscoplasticity. 

B. Time scale defined by the limitations of quasi- 
static analysis. (The upper limit to times considered 
in this class is to some extent arbitrary but could 
probably be limited to a few seconds.) This includes 
times greater than those necessary for the attenua- 
tion and dispersion of elastic waves or the approxi- 
mate solution of problems where wave propagation is 
known to be of minor importance; this class also 
includes most analyses of severe accidents where 
interaction between neutronics and mechanical behav- 
ior is limited to coupling through time derivatives of 
temperature by way of induced thermal displacements 
and the Doppler coefficient. 

C. A broad and somewhat ill-defined time range 
having as its lower end times used in the analysis of 
less severe accidents and as its upper limit a time 
(perhaps minutes) sufficient to include mechanical 
effects noticed during controlled power oscillations, 
reactor shutdown procedures, etc. Perhaps the best 
scale definition would be one equivalent to the time 
associated with low-cycle thermal fatigue. 

D. Very-long-time effects, say from several days 
to years, during which essentially steady-state con- 
ditions have existed. This is the range in which radia- 
tion-induced growth and creep are of greatest relative 
importance. 


Problems of Class A will not be discussed further 
due to their extreme complexity, specialized nature, 
and relative lack of connection to questions of inter- 
est in the normal operation of power reactors. 
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Problems whose time scales fall into Class B are 
almost exclusively related to questions of reactor 
safety and thus fall outside the range of interest in 
this article; one exception would be the use of prob- 
lems of Classes C and D as input to such safety 
analyses. The study of long-time effects as a means 
of specifying the initial deformation and stress field 
in a fuel rod at the time of initiation of a severe 
accident could be used to define a wide class of 
problems whose solutions, although using models and 
methods applicable to Classes C and D, are tailored 
to supply the specific information needed in forming 
a well-posed problem of Class B. 


Discussion 


The two papers!» we will discuss in detail fall into 
Classes C and D. In Ref. 1, time enters only in an 
implicit fashion through various coefficients defining 
radiation-induced creep, swelling due to entrapped 
fission gas, etc.; thus it belongs in Class D. The 
problem treated in Ref. 2 has a range of times 
limited by what could be considered reasonable 
positive reactor periods (i.e., ~30 sec) and thus is 
within the range of Class C. 


MECHANICAL ANALYSIS 


The problem considered in Ref. 1 concerns theo- 
retically predicting the stresses and strains in a 
solid, thinly clad, circular cylindrical fuel element. 
These changes are caused by thermal-cycling growth 
of uranium and its growth, swelling, and creep caused 
by irradiation. Although not explicitly stated, these 
effects are treated without reference to the elastic 
stress state of the fuel and cladding. That is, the 
analysis of all radiation effects is based on purely 
plastic considerations with no reference to, say, 
elastic-plastic boundaries that might exist solely as 
a result of the derived temperature field. As a result 
of this decoupling, elastic material properties such 
as Young’s modulus and Poisson’s ratio do not 
explicitly appear. It should also be noted that, al- 
though practicality dictates that various irradiation 
effects such as those mentioned be treated as sepa- 
rate phenomena, a satisfactory theoretical basis for 
this differentiation would be difficult to outline. Thus 
the separation must be considered to some extent as 
arbitrary. 


Types of Dimensional Changes. As defined in 
Ref. 1, thermal-cycling growth is the gross dimen- 
sional change in the fuel material caused by the 
action of temperature cycles on the crystalline 
structure of the fuel. Such growth is a result of the 
basic anisotropic nature of the crystal structure. 

Irradiation growth, a phenomenon that produces 
Similar gross changes, also results from the anisot- 
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ropy of the fuel but is that portion of the growth due’ 


solely to fuel burnup. 

The term irradiation swelling identifies radiation- 
induced dimensional changes that are caused by 
fission-product gases being trapped within the fuel 
and lowering its density. 

Irradiation creep separates long-term growth ef- 
fects caused by the fuel creep rate being increased 
during fuel exposure to a severe neutron environment. 

Satisfactory physical models for these growth ef- 
fects will, no doubt, eventually become available as 
work in dislocation theory and similar disciplines 
begins to supply methods for relating microscopic 
behavior to observed macroscopic material prop- 
erties. However, at present, one must rely almost 
entirely on empirical descriptions. For example, in 
Refs. 1 and 2, thermal-cycling growth and irradiation 
growth are included by noting that the ratio of final 
to initial length is approximately related to number 
of thermal cycles and percent burnup by an exponen- 
tial relation derived from a fit to experimental data. 
Gross effects due to microstructure have been treated 
in much more exact terms than these in connection 
with the behavior of composite materials (e.g., 
Ref. 8), but such work has had little application in 
the nuclear field. 


Heal-Generation Rate. If the heat-generation rate 
within a particular fuel element is to be considered, 
the local neutron flux must first be specified. If it is 
assumed that the diameter of a single fuel element is 
small when compared with the radial dimension of 
the reactor core, within the limitations of elementary 
diffusion theory, the flux distribution, ¢, in a partic- 
ular element is given by 


g = Aly(81) cos 72 


where /,) = the zero-order modified Bessel function of 
the first kind 
A =a constant for a given element 
B, 7 = constants 
7 = yadial and axial coordinates relative to 
the element center 


~ 
N 


If this flux-distribution equation were used directly 
in an analytical development, the temperature and 
stress fields would be complicated. Therefore, except 
for very local effects, it is usually assumed that 
axial variations in flux are of minor importance. It 
is further assumed that the neutron-diffusion length 
is large relative to the radius of the fuel element, 
thus causing the radial dependence, / (31), to reduce 
to unity. [This latter approximation is quite reason- 
able for aifast reactor system, although it may not be 
true in a thermal reactor. In the thermal system, 
however, the operating conditions are probably not so 
severe, the stresses are likely to be considered less, 
and the resulting errors are probably not of major 














DESIGN AND CONSTRUCTION PRACTICE 15 


importance. The neglect of axial effects allows one to 
consider a plane strain formulation and, hence, 
greatly reduces the mathematical difficulties of the 
problem. Thus the flux has been reduced to a con- 
stant as a first approximation in the analysis, and 
the constant heat-generation rate may be computed 
immediately. After the assumption of constant ther- 
mal conductivity, radial symmetry, and known coolant 
properties at the element surface, the radial varia- 
tion of the temperature can be obtained in a straight- 
forward manner. 


Deformation. In analyses of the radiation-induced 
plastic deformation of the fuel element, it is assumed 
that, in addition to the previously mentioned plane 
strain conditions, the fuel and cladding maintain 
constant densities and that there is no initial or 
plastic-yield-induced anisotropy. Further, both the 
cladding and the fuel are assumed to obey von Mises’ 
yield criterion and to exhibit radiation-induced hard- 
ening that can be described in the same manner as 
conventional strain hardening. With these assump- 
tions the displacement-strain equations, incompress- 
ibility condition, yield criterion, equilibrium, and 
boundary conditions provide sufficient information to 
allow an approximate solution for the state of stress 
in the fuel element if certain empirical material data 
are available. 


Calculation Procedure. The calculation procedure 
may be roughly outlined as follows: Incompress- 
ibility allows the radial displacement to be obtained 
as an integral of the total thermal and radiation 
dilatation. With an empirical description of the 
radiation component of dilatation, which in Ref. 1 is 
assumed to vary exponentially with radial position so 
that the exponential coefficient is a small number 
that allows a simple power-series expansion of this 
term, integration then yields the radial displacement. 
This in turn gives the strain components in the fuel 
and cladding. If the shearing-yield-point stress (i.e., 
the second stress invariant at yield) is known from 
experiment, the radiation-induced strain hardening 
coefficient can be computed in terms of known 
quantities. Then sufficient information is available 
to allow computation of the stress state in the 
cladding and to obtain an iterative approximation to 
the radial variation of the stresses within the fuel. 


Reservations. There are many theoretical reserva- 
tions that could be expressed about the details of the 
computation method used in Ref. 1. However, it is 
perhaps more important to emphasize that all the 
results depend strongly on experimental data on 
effects that may be quite hard to isolate in practice. 
In particular, one should have detailed experimental 
knowledge of the radiation-induced dilatation’ of the 
fuel material and of the behavior of the yield surface 
under irradiation. In addition, even though the radia- 
tion-induced strain hardening of the material is 


specified by a derived quantity insofar as the analysis 
is concerned, it would seem prudent to check such 
derived results against experimental data before 
using this method of analysis as a design tool. 


Evaluation. Overall, because of the numerous as- 
sumptions and the dependence on experimental knowl- 
edge that is difficult to obtain, it seems fair to state 
that serious reservations must be felt about the 
usefulness of these analytical methods to provide 
anything beyond a rough indication of the relative 
Sensitivity of the true state of stress to various 
forms of radiation-induced deformation. 


TRANSIENT GROWTH EFFECTS 

The problem discussed in Ref. 2 is virtually 
identical with that of Ref. 1 with respect to fuel-rod 
geometry and method of analysis. However, the 
object in Ref. 2 is to study transient growth effects 
expected to be of importance during time ranges of 
the order of reasonable positive reactor periods. 
These are introduced into the analysis solely by 
including an exponential time dependence to multiply 
the spatially constant neutron flux. This makes the 
expression for the transient temperature field rather 
cumbersome and, of course, one now works in terms 
of strain rates rather than strains. But, other than 
such obvious modifications, the method of analysis in 
Ref. 1 is essentially repeated in Ref. 2. Hence all the 
assumptions and reservations discussed above re- 
main pertinent. 

In addition, one is faced with the more serious 
problem of justifying the method of analysis in time 
ranges of the order of a minute orless (the numerical 
example in Ref. 2 assumes a reactor period of 36 
sec). It seems unrealistic to attempt to explain radia- 
tion-induced growth effects during such short periods 
without any experimental evidence concerning, for 
example, the plastic flow rules associated with the 
assumed yield criterion. Even if any mention of 
radiation effects is neglected, it is doubtful if the 
equivalent quasi-static plasticity problem for a clad 
cylinder exhibiting strain hardening could be solved 
in the manner outlined. 


Summary 


It seems that an invaluable “ground rule” to apply 
at the beginning of any theoretical transient me- 
chanical analysis of reactor fuel elements is the 
need to include a demonstration that the mathematical 
model chosen applies to the time range under con- 
sideration. Further, unless the only interest of the 
study is to reveal the qualitative relation of various 
hypothesized physical-growth mechanisms to the state 
of stress in the fuel element, the analysis should 
not require experimental material data that are diffi- 
cult or impossible to obtain in practice. 
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OPERATING EXPERIENCE 


Operational Problems and Solutions 


By Rolf C. Skaardal 


Several cases of reactor system component failure 
were reported at the July 1967 Conference on Reactor 
Operating Experience sponsored by the Reactor Op- 
erations Division of the American Nuclear Society 
(ANS). Such failures can be traced to errors in de- 
sign, poor material selection, improper installation, 
and lack of testing; they should be of interest to 
reactor designers and are, therefore, reported here. 
Repair procedures and preparations for repairs as 
well as some equipment modifications are also taken 
into account in this article. 


Failure of Mechanical Components 


Component design and materials problems dis- 
cussed here include the following: 

¢ Hydriding of Zircaloy process tubes in a Hanford 
production reactor 

e Fast-neutron damage in aluminum and beryllium 
core components 

¢ Flow-induced vibration in heat-exchanger tubes 

e Failure and replacement of shaft seals in a pri- 
mary coolant pump 


HYDRIDING OF ZIRCALOY PROCESS TUBES 


Two years after retubing of two of Hanford’s pro- 
duction reactors with Zircaloy-2 process tubes it was 
discovered that they were absorbing significant quan- 
tities of hydrogen.’ The tubes are 1.8 in. in outside 
diameter by 0.037 in. in wall thickness and are about 
47 ft long; cooling-water velocity within the tubes, 
which are surrounded by a mixture of helium and 
nitrogen, varies from 10 to 50 ft/sec. Continued use 
of these tubes appeared in question after one of them 
was removed and found to be very brittle. The ex- 
istence of hydride was confirmed by metallurgical 
and chemical analyses of samples from this tube. 

The hydride layer found on the inside, or water 
side, of the tube was about 0.004 in. thick and showed 
no significant difference in the appearance of longi- 
tudinal and transverse sections. Chemical analyses of 


the entire tube wall showed a total hydrogen content 
in the 1200- to 1400-ppm range, with approximately 
50 ppm in the base metal itself. Recently a tube with 
a total hydrogen content of 4000 ppm was found. 

No apparent direct correlation between total hydro- 
gen content and residence time in the reactor has 
been found. However, the rate at which the hydrogen 
content of the base metal increases seems to bea 
uniform 25 to 30 ppm/year.' 

Once the case layer has been removed, no gradient 
through the tube wall can be detected. The hydride 
platelets, or needles, have never been found to be 
oriented radially except at the Van Stone flange where 
the tube was cold rolled during installation. 


System Conditions. It is common knowledge in the 
nuclear field that hydriding of Zircaloy may result 
from a high-temperature corrosion reaction with 
water or steam. Also, in some cases, it may result 
from a gas-phase reaction where moisture or hydro- 
gen is an impurity in the atmosphere adjacent to 
Zircaloy. But in the Hanford case, the temperature 
of the coolant is too low to be responsible for corro- 
sion reactions sufficient to account for the hydriding 
phenomenon. Also, there is no evidence to indicate 
that the reactor atmosphere could in any way con- 
tribute to the hydride found. 

Laboratory tests have indicated that some hydrogen 
could be introduced into Zircaloy from low-tempera- 
ture mixtures of hydrogen and water vapor.” However, 
the amount of hydrogen which would be introduced 
would be far below the levels encountered here. Ex- 
clusive of the case layer, the maximum hydrogen 
content found uniformly distributed in the base metal 
has been 175 ppm. The dissolving of this much hydro- 
gen normally requires a temperature >300°C. 


Electrolysis, References 2 to 4 indicate that, by 
making zirconium the cathode in an electrolytic cell 
where hot water is the electrolyte, hydrogen can be 
driven into or deposited on zirconium. Thus the 
problem might result from unexpected electric cur- 
rents. However, it did not seem reasonable that the 
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magnitude of the electric currents would be sufficient 
to produce the observed hydriding even though (1) free 
electrons are available from radiolytic decomposition 
of the water, and (2) streaming potentials are caused 
by the high-velocity water. 


Galvanic Coupling. Another possible source of cur- 
rent was the galvanic coupling of the Zircaloy tubes 
with the aluminum-clad fuel elements and aluminum 
dummy spacers. Preliminary tests verified that 
Zircaloy-2 would hydride in hot process water when 
electrically coupled to 6063 aluminum. The prevention 
of current flow, and thus the elimination of this cause 
of hydriding, requires either elimination of the couple 
or insulation of the two metals from each other. 

Dummy elements made of different alternate ma- 
terials, such as nickel-plated aluminum, stainless 
steel, Zircaloy, and carbon steel, were tested in the 
reactor. They all resulted in greatly reduced poten- 
tials or, as in the case of stainless steel, in a couple 
of reversed polarity. 

All these materials, however, have serious draw- 
backs, such as high costs for the first three men- 
tioned, a tendency for severe pitting corrosion of the 
nickel-plated aluminum, and a low corrosion resis- 
tance and high cross section for the carbon steel. 


An Anodized Cure? The practical solution has been 
to insulate the tubes from the aluminum with a high- 
quality anodized coating on the aluminum. Although 
this has been done, detailed operating results are not 
yet available. Anodized coatings will break down 
eventually in hot water, but apparently a good coating 
will not allow significant hydriding during the normal 
residence time in the reactor. Further, anodized 
aluminum has none of the drawbacks mentioned above 
for the alternate metals. 

The oxalic acid anodizing process is used and gives 
a hard scratch-resistant coat of at least 1-mil thick- 
ness. Such a coating has a resistance per square inch 
of approximately 225,000 ohms. After the pieces are 
anodized, they are sealed by autoclaving in water at 
170°C for '/ hr. 


Alternatives. If elimination of the galvanic currents 
fails to control hydrogen buildup, it will be necessary 
to remove the present case layer. Because the pro- 
cess tubes cannot be removed from the reactor, 
several in situ methods have been considered: (1) hot 
vacuum extraction, (2) chemical etching, and (3) me- 
chanical removal by machining, grinding, or grit 
blasting. 

Grit blasting is believed to be the most practical 
means. Trial attempts (both in the laboratory and in 
reactor) to remove given amounts of hydride from 
the tube walls have been fairly successful. Regard- 
less, cleaning an entire reactor would be quite a 
large project; it would take several weeks and use 
about 500 tons of grit. The grit would become highly 
radioactive and could only be used once. 
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NEUTRON DAMAGE TO BERYLLIUM AND ALUMINUM 
CORE COMPONENTS 

Fast-neutron damage and distortion of beryllium 
lattice components have long been recognized as a 
serious problem.',§ Recently, however, similar dam- 
age to Al—Cd shim safety rods (once assumed to be a 
permanent solution to shim safety-rod problems of 
this kind) was found in the Oak Ridge Research (ORR) 
Reactor.’ In addition, two aluminum lattice pieces 
that held experimental facilities were also found to 
be distorted. 


Bowing. In each case, bowing occurred— appar- 
ently from irradiation damage— with effects like 
those observed for fast-neutron damage to beryllium. 
This assumption is strengthened by test results of 
embrittled aluminum cladding from irradiated HFIR 
target rods, several of which have failed due to cracks 
in the cladding. In each instance at the ORR, the 
dimensional changes were sufficient to cause either 
improper fit of the component or interference with 
adjacent lattice pieces, or both. 

The Al—Cd shim safety rods, located in lattice 
positions F-4 and F-6 of the ORR lattice (Fig. 1), 
failed to pass a manual check for freedom of move- 
ment prior to the reactor startup on Sept. 20, 1966. 
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Fig. 1 The Al—Cd shim safety rods in positions F-4 and 


F-6 of the ORR lattice failed a manual check and subse- 
quently were found to bow toward the fuel region.’ 
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Figure 2 shows the normal alignment of these rods 
during operation. When the rods were checked with a 
straightedge they proved to be bowed toward the 
largest source of neutrons (the fuel region) as were 
Be—Cd rods that had been in the same lattice posi- 
tions earlier. 
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Fig. 2. Cross section of the ORR showing the alignment of 
the F-4 and F-6 Al—Cd shim rods.’ 


Veasurements. Distortion of the aluminum portion 
of the rods (which were sectioned) was measured with 
a contour-measuring apparatus. This apparatus, Fig. 
3, consists of a calibrated slide bar, a machined 
structural bed plate, a dial indicator with mounting 
clamps, and two vises for positioning the parts to be 
measured. Calibration curves were obtained with bar 
stock known to be flat to 0.0005 in. Measurements of 
deflection were noted at 1l-in. intervals parallel to 
the long axis of the aluminum section, and the data 
were adjusted from calibration curves. 

Figure 4 shows the east-west measurements for 
the F-4 rod plotted as a function of distance from the 
top of the fuel region. (North— south deflections were 
relatively low, as expected.) The slight bowing to- 
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Fig. 3 Device used tomeasure bowing in ORR components.* 


ward the north is expected for the rodin position F-4, 
and similar results were found for the rod in F-6 
because fuel was located in position F-7 during part 
of the irradiation period. 

It is interesting to note that, for measurements on 
opposite sides of each rod, the maximum deflection is 
5 to 10 mils greater on the side nearer the source. 
This characteristic has also been exhibited by dam- 
aged beryllium. Growth would be greatest in the 
region of highest flux and could indeed cause the 
aluminum (or beryllium) sections to be thicker near 
the center than at the ends. 
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Fig. 4 Bowing measurements made on two sides of the 
Al—Cd rod in core position F-4.° 


Table 1 lists observed deflections and estimated 
fast-neutron fluences. The data show that aluminum 
in core positions F-4 and F-6 bows at a rate of about 
0.015 in./year—the rate for beryllium in the same 
positions is ~0.010 in./year. 

Table 1 also contains measurements of bowing in 
two aluminum lattice pieces. These items were lo- 
cated on the periphery of the fuel region, specifically 
in lattice positions A-7 and D-8. One is essentially a 
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Table 1 IRRADIATION HISTORIES AND BOWING 
MEASUREMENTS? OF ALUMINUM LATTICE 
COMPONENTS IN THE ORR* 


Fluence,tf 


Lattice >1.0 Mev, Bowing ,t 
Component position Mwd neutrons/cm? in, 
Al—Cd shim rod, 
Al-3S F-6 21,491 9.97 x 10?! 0.039 
Al—Cd shim rod, re 
Al-4S F-4 21,491 9.22 x 10° 0.030 
Hydraulic tube No, 2 D-8 28,711 1.1 = 1074 0.029 
Np tray A-7  ~67,549 1.8 x 1072 0.060 


*The ORR produces ~9400 Mwd(t)/year and operates at a nom- 
inal power of 30 Mw about 86% of the time. 

iEstimates are based on measurements at the outer edge of a 
fuel element adjacent to a beryllium reflector piece. 

tBowing measurements given are the maximum measured 
values that occur on the side adjacent to the fuel region. 


hollow aluminum box, and the other is an aiuminum 
block with several drilled holes. 

The hollow aluminum box (the N, tray) was bowed 
in two directions as expected. The side away from the 
fuel bows less than does the side adjacent to the fuel, 
again indicating a “fattening” of the middle section of 
the lattice piece. 

Damage of beryllium lattice pieces due to irradia- 
tion was established some 5 years ago in the ORR.* 
(A similar case of beryllium distortion occurred at 
the Materials Testing Reactor (MTR) where the 
beryllium section of a shim rod was found to be con- 
siderably bowed after 27,193 Mwd of irradiation. °) 
Table 2 gives the irradiation history and bowing mea- 
surements of damaged beryllium lattice components 
in the ORR. 


fable 2. IRRADIATION HISTORIES AND BOWING 
MEASURE MENTS® OF BERYLLIUM LATTICE 
COMPONENTS IN THE ORR 


Lattice 


Fluence,* 





posi- >1.0 Mev, Bowing,f 
Component tion Mwd__ neutrons/cm? in. 
Be—Cd shim rod, Be-1S F-6 24,323 1022 0.031 
Be-—Cd shim rod, Be-2S  F-4 25,432 ~— 10 0.028 
Be reflector, R-10 F-5 36,975 1.8 x 1072 0.056 








*Estimates are based on measurements at the outer edge of ¢ 
fuel element adjacent to a beryllium reflector piece.” 

*+Bowing measurements given are the maximum measured values 
that occur on the side adjacent to the fuel region. 


HEAT-EXCHANGER TUBE VIBRATIONS 


Flow-induced vibrations and resulting damage of 
heat-exchanger tubes have in the last few years been 
recognized as a serious problem in need of thorough 
engineering attention. The most well-known case in 
this country probably concerns the Advanced Test 
Reactor (ATR).! Similar High Flux Beam Reactor 
(HFBR) experience has been reported.® However, in 
the HFBR case, the symptoms were recognized and 
the problem solved before any tube rupture occurred. 


HF BR-Exchanger Design. The two HFBR primary 
coolers (Table 3 and Fig. 5) are U-tube counterflow 
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exchangers with the heavy water on the tube side and’ 
cooling water on the shell side. Rattling noises were 
heard in these large coolers (68.5 x 10° Btu/hr) dur- 
ing an early test at rated secondary flow. Design 
Similarities in size and baffle spacing which existed 
between these coolers and those of the ATR precipi- 
tated a study to determine the effect of reduced 
cooling-water flow. It was found that at about one-half 
rated flow the rattling disappeared. 

Analysis indicated that the natural vibrational fre- 
quency of the tubes was sufficiently close to the 
driving frequency at rated cooling-water flow to 
cause some resonant vibrations. Because of the 
relatively long unsupported tube spans, the U bends 
were particularly suspect. Therefore it was decided 
to undertake a vibration survey to define the problem. 


Table 3.» DESIGN DATA FOR THE HFBR PRIMARY COOLERS* 





Tube side Shell side 





Fluid D,O Cooling-tower 
water 
Velocity, ft/sec 7.3 4.65 
Rated flow, gal/min 9000 4000 
Pressures 
Operating, psig 250 12 
Drop at rated flow, psi 10 12 
Temperatures at rated flow 
D,O in, °F 134 
D,O out, °F 120 
Cooling water in, °F 85 
Cooling water out, °F 119 
Number of U tubes 1350 
Tube specifications ¥/, in. in outside diameter 
x 16 BWG 
Tube-to-tube clearance, in. 6 
Tube pitch, in. Al), 
Length between baffles, in. 29//, 








D20 outlet 
Cooling - water outlet 
D20 inlet 
Fig. 5 Vibration in the as-designed HFBR primary cool- 


ers proved excessive; peak-to-peak amplitudes reached 
0.150 in. at vated 4000 gal/min flow.® 


OPERATING 


Vibration Survey. Vibration-sensing elements 
were inserted into the previously drained and dried 
tubes of the cooler. The sensing elements were precal- 
ibrated piezoelectric-type accelerometers that could 
sense vibrations in two mutually perpendicular direc- 
tions. The signals from the accelerometers were 
amplified, integrated, and recorded on an oscillograph. 
Twenty-nine tubes in one quadrant of one cooler were 
examined; measurements were made at several posi- 
tions in each tube. Some of the tubes were selected in 
a random fashion— others because of the greater 
likelihood of vibrations (e.g., tubes near the cooling- 
water inlet nozzle). Longitudinal positions selected 
were uSually the midpoints between the baffles. Mea- 
surements were not made in the curved section of the 
U tube because of the inaccessibility. 

In tube positions where the unsupported length was 
two baffle spacings (~ 60 in.), the vibrational fre- 
quencies were 30 to 40 cps. Between the last baffle 
and the curved section of the U tube, however, the 
frequencies observed were between 5 and 15 cps. 
Where unsupported length was ~30 in., frequencies 
approached 140 cps. 

View ports were cut into the side and U-bend end 
of the shells. Vibration visible through side ports did 
not appear serious. However, through the end ports, 
three of four tubes could be seen vibrating with sub- 
stantial amplitudes. Motion pictures indicated ampli- 
tudes as high as 0.15 in. peak to peak at a flow of 
4000 gal/min. No vibration was noted at 2000 gal/min. 

The conclusions were that (1) there was little 
likelihood of dangerous vibrations in the straight sec- 
tions of the tubes, (2) there was serious vibration in 
the U bends at rated flow, and (3) tube vibration at 
one-half rated flow appeared insignificant. 


Revised Design. Reduction of tube vibration in the 
coolers without seriously limiting the heat-removal 
capability (it was necessary to increase primary-to- 
secondary log mean temperature difference by ~25%) 
was accomplished by cutting new outlet holes in the 
U-bend end of the coolers and splitting the cooling 
flow so that half the rated flow went through each 
shell pass. The goals were to restore cooling flow to 
its original value, to limit shell-side velocities to 
one-half of the original value over the straight section 
of the tubes, and to almost eliminate flow over the 
U bends (Fig. 6). 

The modifications cost about $29,000 (20% of the 
initial cost of the exchangers) and were successful; 
at rated flow, no tube vibration has been detected. 

Tube damage due to flow-induced vibrations in heat 
exchangers is a serious problem which has been given 
too little design consideration up to now. The HFBR 
case is another strong reminder that this must be 
done. Improved designs, analytical approaches, and 
in-shop testing procedures are being developed. For 
the new Argonne Advanced Research Reactor facility 
(for example), three 30-Mw heat exchangers are at 
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D20 outlet 
Cooling — water inlet 
D»0 inlet 
Fig. 6 Modified primary coolers for the HF BR almost 


eliminate flow over the U bends.§ 


present under contract negotiations. A full-scale 
mockup test is specified, and this contract will stipu- 
late a manufacturer’s guarantee against flow-induced 
vibrations. 


CIRCULATION -PUMP SEAL FAILURE 


When the Big Rock Point Nuclear Plant was re- 
started after a routine shutdown in November 1966, 
it became apparent that the shaft seals in the No. 2 
circulation pump had failed.® In 1964 these controlled- 
leakage, mechanical shaft seals were changed from 
the original rotating hard-face Stellite to aluminum 
oxide (sapphire) for the No. 1 pump and to tungsten 
carbide for the No. 2 pump and had accumulated 
~16,000 hr of operation when the tungsten carbide 
seals failed. 


Planned Repairs. A schedule was laid out for re- 
placing the failed seals in situ, with estimates of the 
time and number of personnel required to accomplish 
each step. On the basis of radiation levels (up to 2 r/ 
hr in the working area) known from surveys conducted 
during the refueling outage, it was estimated that the 
job would require 16 men and could be done in about 
12 hr of working time. 

Repairmen were imported, briefed on radiation 
protection practices at Big Rock Point, and assigned 
a specific step in the seal-replacement procedure. 
They then studied their specific jobs with site per- 
sonnel. Replacement parts were brought out, and, 
with drawings and pictures, each job was described 
in detail. 

The work started at 8:00 a.m. on a Saturday, and 
the plant was expected to be back in full service by 
Monday. No problem was anticipated— contingencies 
had been included. But... 
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Complications. When the outer seal unit was ex- 
amined, it was found that the tungsten carbide rotating 
seal face had shattered and caused severe damage to 
the pump shat and other parts inthe seal mechanism. 
This meant that the shaft would have to be removed 
and repaired or replaced. The magnitude of this job 
dictated that the entire pump be removed from the 
high-radiation pump room. The shattering of this seal 
is believed to have been caused by an insufficient 
seal-support arrangement, which later was redesigned 
by the pump manufacturer. 

Radiation-Exposure Problems, The work required 
to remove the pump from the recirculating pump 
room was straightforward and did not demand a de- 
tailed briefing. When the pump was lifted free of its 
casing, the radiation reading was 20 r/hr at the 
impeller. The original location chosen for storing the 
pump was not suitable for this high radiation level; 
therefore the pump had to hang on the reactor crane 
until a suitable location was readied. While a pit on 
the refueling deck was being prepared, long-range 
decontamination techniques were employed to reduce 
the radiation reading at the impeller; this was low- 
ered to approximately 9 r/hr but no further. The 
pump was then placed in the pit. Sheet lead was 
used to reduce the radiation field in the area. The 
men who had worked on this phase of the job had ac- 
cumulated individual exposures of 1.8 to 2.7 r. 

A stand was designed and fabricated to hold the 
pump in its normal upright position and also in an 
inverted position for impeller removal. The pump 
and stand were submerged in water (in a site- 
fabricated steel tank) to obtain a radiation shield and 
prevent airborne contamination. Sheet lead was also 
used where practicable during the pump disassembly 
to reduce the radiation field. 

The job of disassembly and assembly required 
3 weeks and was done by three groups of three men 
each. Each group worked 1 week on the job and ac- 
cumulated an average exposure of 2.5 r in 40 hr of 
work. A plant engineer who supervised accumulated 
approximately 1.5 r during the 3-week period. 


Pump Improvements. During reinstallation of the 
pump, several areas were noted in which modifica- 
tions could save considerable time and exposure: 

e Swagelok or similar fittings should be used on 
the seal-system controlled-leakage lines. 

¢ Connectors should be used on the instrumentation 
leads. 

¢ Provisions should be made to simplify draining 
and refilling bearing-oil reservoirs. 

e Appendages should be minimized to ease lifting 
and lowering through tight spaces. 


Work in Radiation Areas. Maintenance workon 
highly radioactive equipment or on equipment in high- 
radiation areas at Big Rock Point and other facilities 
leads to the following general conclusions: 
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(1) Procedures must be flexible but as detailed as. 
possible. 

(2) Skilled manpower with the required radiation- 
exposure history must be available. (As more of the 
utility's generating capacity becomes nuclear, the 
pool of outside personnel with low radiation exposures 
will decrease.) 

(3) Equipment should be designed to facilitate 
maintenance. Plant layout must be aimed to reduce 
needless exposures that come from poor and im- 
properly placed equipment. 

(4) Key personnel should conduct a follow-up study 
immediately following each major job. Problem areas, 
solutions, and suggested solutions should be recorded. 


Reactor Modifications 


Ultimately, difficulties with components and mate- 
rials lead to major reactor modifications. 


PREPARATIONS FOR REPLACING 
A REACTOR VESSEL 


Since first criticality in 1957, the NRU, a heavy- 
water-cooled and -moderated reactor (Fig. 7), has 
experienced significant corrosion in the aluminum 
reactor vessel.'? As a result, much of the vessel— 
which was site assembled— may have to be replaced. 
Detailed plans have been developed and trial runs of 
removal and disposal operations have been made, 
often with full-scale mockups. A PERT program has 
been developed to make maximum use of manpower 
and resources. The program calls for full-time ef- 
fort— 24 hr/day, 7 days/week— during a reactor 
downtime of about 6 months. The prognosis, of 
course, is for successful completion. 


Corrosion Damage. The NRU reactor-vessel alu- 
minum sidewall is a rather complicated structure; it 
includes an annulus immediately outside the heavy- 
water vessel and, outside that, a light- water reflector. 
Corrosion is caused by water that leaks into the an- 
nulus from the reflector and from water-cooled as- 
semblies irradiated in the annulus. An 8-ft-square 
aluminum window Y,, in. thick, used in lieu of the 
reflector in front of the graphite thermal column, has 
completely corroded through. If the sidewall of the 
reactor vessel should fail, attempts will be made to 
repair it in situ. Certain repairs, however, will not 
be possible. 


Vessel Sidewall Delails. The cylindrical aluminum 
sidewall, which is 0.188 to 0.313 in. thick, is sealed 
between top and bottom headers of stainless steel. 
The sidewall is secured to the two header flanges by 
72 bolts at the top and 72 bolts at the bottom. These 
bolts shoulder on the headers and fasten to inserts 
in the sidewall flanges. To form a double ‘seal, a 
metal-to-metal butt seal and a stainless-steel O- 
ring are used. 
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Fig. 7 Because of corrosion the aluminum reactor vessel in the Canadian NRU reactor may be replaced; 


preparations have been made.'" 
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Another aluminum sidewall assembly has been de- 
livered. Although similar to the original, it incorpo- 
rates a new seal section made up of two '/p-in.-wide 
aluminum rings, one on each side of the bolt-circle 
centerline. These rings are integral parts of the 
aluminum flange and will mate to the seal surface of 
the existing stainless-steel header. Part of the pre- 
paratory program has been to ready this assembly. 

Because of high-radiation fields, all operations 
must be conducted from behind heavy shielding. 
Planned facility design makes it possible to remove 
and retain behind shielding such large units as the 
bottom shield and bottom header. These units are 
anticipated to have contact radiation levels in excess 
of 20,000 r/hr. 


Removal. Figure 8 shows the basic steps in the 
procedure for removing the sidewall assembly from 
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Fig. 8 Steps in the removal of the sidewall and lower 
header” 


the reactor. The 144 bolts at top and bottom must be 
removed. The upper shields and header remain in 
place as part of the reactor structure. Next, the 
aluminum sidewall and the bottom header and shield 
are again bolted together with three long capstan 
bolts that are operated from above. Then the as- 
sembly is (to be) lowered to a carriage. 

Prior to these operations a steel skirt with a solid 
top equipped with a lifting yoke and fitted with an 
open-ended nylon bag will be hung at the ceiling of 
this room. When the carriage is properly positioned, 
the skirt will be lowered over the aluminum sidewall 
assembly and the lifting adapter secured tothe vessel. 
After the floor slabs are lifted away, the disposable 
portion of the vessel will be enclosed in a metal and 
nylon container (for contamination control) and this 
package will be removed from the shielded area to 
the disposal skid positioned in the reactor hall. Op- 
erations will be viewed on closed-circuit television. 
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Disposal, One hundred days after reactor shut- 
down, the radiation fields from this 29,000-lb package 
are estimated to be about 350 r/hr at 25 ft—the 
distance behind a shielded vehicle at which it will be 
towed to a disposal area. The preferred vehicle to 
tow the package is a mobile crane equipped to hitch 
and disconnect remotely by hydraulic means. The 
disposal area has been prepared so that the tow 
vehicle can drive in, disconnect, and drive out, leav- 
ing the skid (a solid-tired flat-bed trailer) in the 
disposal pit. 

Trial runs of the disposal procedure have been 
completed with a dummy vessel under winter and 
summer conditions. 


Custom Tooling. Much special tooling is required. 
The most interesting is a very complex $100,000 
surface-conditioning tool. Because the condition of 
the sealing surfaces on the stainless-steel headers is 
unknown, the conditioning tool must be able to in- 
spect, brush, hone, grind, and clean as required to 
effect a 16-,in. surface finish. Now being tested, the 
tool is a three-armed spider with a vertical main 
drive stem. Tooling and inspection heads for various 
functions attach to each arm. A borescope viewing 
head and vacuum-cleaning facilities are provided. 

The surface-conditioning tool is designed to oper- 
ate in an environment of about 20,000 r/hr through a 
4.7-in.-diameter vertical hole and to work on sur- 
faces about 11 ft in diameter, 28 ft below the manned 
level. 


INSTALLATION AND TESTING 
OF A STEEL POOL LINER 

The Texas A&M University swimming-pool reactor 
facility has had pool leakage problems since the 
initial construction in 1961 (Ref. 11). Failure of the 
epoxy coating in the main pool section and the alu- 
minum liner in the stall section was the result of 
improper installation and lack of adequate testing pro- 
cedures. A solution was sought by contacting other 
pool-reactor operators concerning their liner experi- 
ence.'* When the pool leakage flow reached 1000 gal/ 
day in March 1966, the decision was made to install 
a welded stainless-steel pool liner. Installation and 
related modifications were completed in a 4-month 
period at a cost of about $51,000; design cost is esti- 
mated at an additional $5000. (Some costs for pool 
and piping improvements are not included.) 


Design. The walls of the 304 stainless-steel liner 
are constructed of 10-gauge sheet, and the floor is of 
/,-in, plate. Drain channels in the concrete beneath 
the floor of the liner collect any leakage. Water from 
these channels drains into the radioactive waste- 
disposal system. 

“Hat channels,” Fig. 9, were formed from 10-gauge 
material and positioned vertically on the walls to 
facilitate vertical welding of the wall sheets. These 
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Fig. 9 “‘Hal channels,” analogous to furring strips, were 
installed vertically on the concrete walls of a reactor pool 
to facilitate installation of a steel liner"! 


channels are held in place by locknuts on powder- 
actuated threaded studs. Vertical slots in the hat 
channels accommodate the studs and allow slight 
movement along the studs. Strips of 10-gauge mate- 
rial were recessed in the pool walls as backup plates 
for the horizontal welds. The base molding of the 
liner is formed by an angle */,, by 2'/, by 2!/, in., and 
this molding serves as an expansion joint between the 
floor and wall sections, Fig. 10. A similar angle 
trims the top of the liner. 


° 


yaar 
\ 
\ 
| 


Ane 
a AEIe 





N 
N 
iN 


arn aes aN SR 7 


V/A 





i TS ial ES I. LES 























WY . : 
//} IS, 7 zz 
Yi Ma. 


Fig. 10 
formed by welding a*/g-in.-thick angle to both surfaces. 





An expansion joint between pool wall and floor is 
: 11 


Installation. The floor was the first section in- 
stalled. Sheets of '/;-in. plate were positioned and 
welded. Vertical hat channels were installed over the 
old epoxy coating. The 10-gauge preformed wall 
sheets were positioned and tack-welded to:the hat 
channels. Gussets were fitted and welded to steel 
penetrations and to the liner. For cost and strength 
reasons, the aluminum gate jamb between the two 


sections of the pool was not replaced. Rather, stain- 
less-steel sheets were formed around the gate and 
welded to adjacent wall sheets. 


Welding. Final welding was done with a short-arc, 
automatic-feed, inert-gas welding machine after all 
wall sheets had been tack welded in position. This 
method proved to be much faster than rod welding; 
however, it had the disadvantage that the amount of 
weld material deposited was sometimes inadequate; 
when the welds were ground flush (to improve the ap- 
pearance of the liner) some welds were too thin. The 
short-arc machine was chosen to limit the heat 
applied and thus reduce buckling. Regardless, con- 
siderable buckling of the hat channels occurred at the 
top and bottom of the liner, but this buckling was cor- 
rected (later) by installation of the base angle and top 
trim angle. Floor plates were rod welded. 

A rod-welding technique should be used throughout 
because the buildup is better and much less splatter 
is produced. Both factors significantly reduce the 
amount of grinding required and produce more satis- 
factory results. Where grinding is not required, the 
rod weld also produces the better finish. 


Testing. First testing of the welds was done with 
a vacuum bubble-box method—a technique that was 
not successful. Successful testing was accomplished 
with a Freon-22 leak-detector system. Gas was 
introduced through the pool-leakage drain system. 
After large leaks around pool penetrations were re- 
paired, several smaller leaks were discovered. After 
several inspection and repair cycles, a raw-water 
test indicated a leak rate of about 25 gal/day. 


Performance. The performance of the liner has 
been excellent. In addition, facility appearance has 
been improved. Experience has confirmed that pool 
leakage is only 25 gal/day. Since the evaporation rate 
is on the order of 200 gal/day, no attempt has been 
made to correct the leakage; it presents no opera- 
tional problem. 
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FLUID AND THERMAL TECHNOLOGY 


Thermophysical Properties 
of Sodium—Recommended Values 


By Gerald H. Golden, John V. Tokar, and David Miller 


Comparisons of safety studies of sodium-cooled fast 
breeder reactors will be more meaningful if the 
various groups making such studies use the same 
values for the thermophysical properties of sodium. 
So that “best values” of these properties could be 
recommended, the literature on the thermophysical 
properties of the liquid and vapor phases of sodium 
has been evaluated.' The resulting recommended 
values are reported here. 

Interest in sodium as a heat-transfer fluid and as a 
working fluid for power cycles has triggered numer- 
ous studies of its properties. In turn, critical com- 
pilations have been necessary every few years, ?~"" 
The newest compilation, which attempts to update 
through 1966, includes: 


Liquid density 

Vapor pressure 

Thermodynamic properties (including P-V—T) 
Transport properties (conductivity, viscosity) 
Electrical resistivity 

Surface tension 

Liquid thermal-expansion coefficient 

Liquid compressibility and sonic velocity 


In each case the results of different investigators 
have been compared, and recommended “best” values 
are presented. The recommended values are pre- 
sented in tabular form to facilitate engineering use 
of the property data. Unabridged versions of the 
thermodynamic charts and tables are available in 
Ref, 1. 

The temperature, pressure, and specific volume at 
the critical point are not included in the compilation 
because of the considerable uncertainty about the 
values reported. However, recent estimates of the 
critical constants and equations for application in the 
critical region are given by Miller’! for sodium and 
the other alkali metals. 


Density of Liquid Sodium 


Table 1 lists the methods that various investigators 
used to measure the density of liquid sodium, the 
temperature range over which they measured it, and 
the average extent by which each set of data differs 
from values calculated by the equation 


p,(Ib/cu ft) = 59.566 — 7.9504 x 10-* 
- 0.2872 x 107°? + 0.06035 x 10°73 = (1) 
where the temperature, /, is between 208 and 2500°F. 


Equation 1 was derived by Stone!” on the basis of the 
data of the investigators listed in Table 1. 


Table 1 SUMMARY OF DENSITY MEASUREMENTS 
FOR LIQUID SODIUM 








Temp. Av. deviation 

Investigators Method range, °F with Eq. 1, % 
Ewing! Dilatometric M.p. to 503 + 0.08 
Jackson*? Buoyancy 937 to 1314 —0.74 
Rinck*® Buoyancy 804 to 1183 +0.15 
Hagen*® Dilatometric M.p. to 336 + 0.05 
Stone!? Dilatometric 1577 to 2491 +0.17 

(pycnometers) 
Novikov Buoyancy of 

Steel sinker 248 to 505 —0.14 

Tungsten sinker 275 to 1324 +0.71 
Nishibayashi®*! Buoyancy 486 to 1580 + 1.09 





Corroboration. The values obtained with Eq. 1 
were corroborated by two other investigations. An 
equation derived by Thomson and Garelis® for the 
208.1 to 1184°F range agrees with Eq. 1 within 0.1%. 
Goltsova!? used a molybdenum pycnometer for mea- 
surements at 212 to 2640°F, from which she derived 
an equation that agrees with her data within an 
average of 0.5%; her equation agrees very well with 
Eq. 1 over the range 212 to 2500°F (at 2500°F, den- 
sity from the two equations agrees within 0.39%). 
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Figure 1 plots Eq. 1 and the results of two other 
investigations. *:*4 The equations derived by Tepper" 
and by Weatherford® do not agree as well with Eq. 1 
as do the equations derived by Thomson and Garelis® 
and Goltsova.'® As the source of the graphical data he 
uses, Weatherford® cites Meisl and Shapiro, ° who in 
turn cite Miller,’ who cites Ewing’® as his source. 
However, in the experiments, Ewing’® measured 
sodium density by a dilatometric technique over the 
range 225 to 503°F; apparently, the other authors’ ’® 
extrapolated the original 225 to 503°F data to 1292°F 
and 2240°F. In view of this large extrapolation, and 
on the basis of the close agreement of the work of 
Thomson and Garelis,*® Goltsova,'* and Stone, Eq. 1 
is recommended as the best representation of the 
density of saturated liquid sodium over the range 
212 to 2500°F, 
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ture. The curve of Stone’* and Goltsova'’ is recommended. 


Saturation Vapor Pressure 


Comparison of various vapor-pressure equations 
for sodium is complicated because some equations 
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are fits to specific sets of experimental data, whereas 
others are “pooled” results. 


FITS TO DATA 


Stone’? measured the saturation vapor pressure of 
sodium with a static-capsule technique. A least- 
squares fit of the data taken in the range 1624 to 
2539°F resulted in the Kirchoff-type equation 


9,980.94 
T(R) 
— 0.61344 log 7(°R) = (2) 


log P,(atm) = 6.83770 — 


This equation gives anormal boiling point of 1618.6°F. 
The experimental data deviate from values computed 
with Eq. 2 by an average of 0.37%. An independent 
“third law check” demonstrated the internal con- 
sistency of the vapor-pressure data. 

Table 2 summarizes the results of related inves- 
tigations since 1955. Three of the six equations give 
normal boiling points that agree to within 0.5°F of 
that given by Eq. 2. Figure 2 compares these data 
and shows that four of the six equations agree with 
Eq. 2 to within 6% in the range 1270 to 2530°F. The 
Makansi'® equation agrees best. 


POOLED CORRELATIONS 


Table 3 summarizes the ranges of applicability of 
the vapor-pressure correlations, the equations, the 
predicted normal boiling point, and lists the refer- 
ences cited by each investigator. Figure 3 plots the 
deviations of six “pooled equation” results from the 
values predicted by Eq. 2; extrapolations are over 
wider temperature ranges than intended by some of 
the authors. As shown in Fig. 3, below 1600°F, most 
of the pooled correlations differ rather uniformly 
from Eq. 2. On this basis (and considering the close 
agreement between the equations of Stone’? and 
Makansi'® above 1600°F), we recommend the Ditch- 
burn and Gilmour’ equation below 1600°F and the 
Stone!” equation for higher temperatures. A similar 
choice was made earlier by Makansi, Selke, and 
Bonilla. '® 


Thermodynamic and P-V-T Properties 


Sodium vapor is composed of the monomer, dimer, 
and, reportedly, tetramer; the fraction of each spe- 
cies depends on both temperature and pressure. !” 
Saturated vapor contains about 10 mole % dimer at 
1500°F and 18 mole % dimer at 2500° F. If the ideal 
monomeric gas is taken as the base, any deviations 
from ideality in the vapor are caused bythe effect the 
higher molecular-weight species has on the mean 
molecular weight of the vapor and on the nonideal 
compressibility and mixing effects of the individual 
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Table 2 VAPOR-PRESSURE INVESTIGATIONS SINCE 1955 
Temp. Derived equation Derived normal 
Method range, °F for log PB, atm boiling point, °F 


Investigators 





Bowles and Rosenblum** Static capsule 


Sowa?? Evaporation from crucible 
at fixed pressure 


Kirillov and Grachev*! Static capsule 


Makansi!* Refluxing capsule 


Achener and Jouthas®® Flash boiler 


Fischer’? Quasistatic variation 


Elastic-membrane- sealed 
static capsule 


Vinogradov and Volyak* 


2 9,363.82 
7 3417 961. — = a 31: 
1470 to 3417 4.51961 TOR) 1612.1 
Cc Cc a+ J 
1651 to 2533 1.56719 — 2:491.55 1618.5 
1616 to 2372 Authors’ smoothed datat 
Q9290KEf 
1148 to 2075 1.521 1618.6 
Po 9,770.25 
1129 to 1752 4:71924 — 2110.29 1610.6 
T 
971 to 1484 4.6279 — 2:616.14 1618.2 
T 
1096 to 1628 Authors’ smoothed data § 1621.1 





*Derived from least-squares fit to data. 

'Table 1 in Ref. 54. 

tFit of equation of Kirchhoff-type tried and found to be 
§Table 3 in Ref. 56. 


species. Before the P—V—T measurements became 
available for sodium vapor, its mass specific volume 
at moderate pressures was computed from ideal-gas 
thermodynamic properties, generally neglecting the 
compressibility and mixing effects. 


HEAT OF DIMERIZATION 


An outline of the thermodynamics of sodium pro- 
vides a basis for discussing techniques that have 
been used to analyze pertinent experimental data. 
Initial assumptions are that the vapor is saturated, 
consists only of the monomer and dimer, and that 
these are ideal gases. For the dimerization reaction 
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(see Table 2) from those recommended by Stone.!* Portions 
of curves based on extrapolations are shown by brokenlines. 


no better than two-term fit (Ref. 18). 


2Na(gas) = Na,(gas) 








AF} = RT In ky (3) 
where 
a T r 
Re ii ma = No = a... am (4) 
P2, P,Ni P,(1—N,)° 
Dividing Eq. 3 by 7 and rearranging 
-R ln ko = —s = (£4) 
T T Je 
-2( i) AH, (5) 
T 1 = 


For convenience, let (AHj))»>= AD§, the heat of di- 
merization of the ideal gas at absolute zero. 

The thermodynamic properties of ideal gaseous Na 
and Na, were computed by Benton and Inatomi!® and 
later by Evans.”’ The values of the free-energy 
functions [(F°~—H))/T], and [(F°—H§)/T]2 reported 
by these two groups are in excellent agreement and 
in accordance with the tables of Stull and Sinke”! 
and with the JANAF tables.” 

The value of AD) has been determined experi- 
mentally by a number of investigators, as described 
in detail by Thomson and Garelis.° The values pre- 
sented range from —8400 + 100 to —9300 + 500 cal/ 
gram atom Na. The vapor composition is extremely 
sensitive to the value of AD5. 


HEAT OF VAPORIZATION 


For vaporization to the ideal monomer 


Na(liquid) = Na(gas) 


AF° =—RT In Py, (6) 


z 
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Table 3 VAPOR-PRESSURE POOLED CORRELATIONS 








Normal 
Temp. Derived equation for boiling 
Investigators range, °F log P,, atm point, °F Refs. 
; as 17 - : 10020.6 " ja se ats _— 
Ditchburn and Gilmour'' M.p. to 1790 6.4818 — os 0.5 log T 1617.5 57-68 
a F eee OTS ; Bs aia i 
Bowles and Rosenblum” 1149 to 3418 4.54025 — a 1618.4 16, 53, 69 
5 9448.38 
T 3 Jarelis® 26 25 Ace. .) oo 
Thomson and Garelis 260 to 1250 4.57791 To 14.65 
» foe he 8718.3 ie ig — 
1250 to 2240 4.36383 — ~=—sa 7 1621.3 16, 57, 59, 60, 
T — 83.16 . 
63, 64, 66 
: 0050 
Meisl and Shapiro® 240 to 1040 8.286 — 10050 — 0.698 log T 
10104 
Weatherford?® 1040 to 2240 8.402 — - - 0.723 log T 1630.3 
10 _ ‘ 10395 ee tea a es dle 
Kelley’ 357 to 1620 9.1882 — ——— — 1.274 log 7 1637.3 97-66, 71-75 
9 9980.9¢ 4 4 
Stone! 1624 to 2539 6.8377 — — ~ 0.61344 log Tt 1618.6 (Eq. 2) 
Hicks” M.p. to 2070 See Ref. 54 1620.9 16, 57, 58, 63 





*Recommended equation for vapor pressure at =1600°F. 
‘Recommended equation above 1600°F. 
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below 1600°F - . 12 0 1600°F 
€ j F and those of Stone’ above 16 F. for enthalpy of liquid sodium.”4 


Thus in Fig. 4 is a curve derived from an equation by 
, To — pe ae? Shpilrain™* for the enthalpy of liquid sodium. Dis- 

(AH®), = -RTIn Py, — T\(4—— 9) - (= 40)) (7 casera Siig be : 
Na T : crepancies in the sodium vapor-pressure data stimu- 


i T 





lated a number of analyses‘:'®.”° to derive values for 

where (AH;),, is the heat of vaporization to the ideal A Dy; and (AH)... 
monomer at absolute zero. The heats of vaporization as computed by the 

Ginnings, Douglas, and Ball”? measured precisely virial method with Eq. 8 and by the quasichemical 
the heat capacity of condensed-phase sodium in the method with Eq. 9 are shown in Fig. 5. 
range 547 to 2106°R, from which they derived em- 
pirical equations for the enthalpy, specific heat, and aie = Kr; ea 7 aPs (8) 
entropy of condensed-phase sodium from 492 to p p,) aT 
2111°R. Stone! measured specific heat at various 
points between 1062 and 2610°R and extrapolated the se = Et ee ithe (9) 
Ginnings correlations to 2610°R (Fig. 4). Also shown ; M 
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Fig. 5 Heat of vaporization of sodium as computed with 


virial and quasicheniical methods. 


Figure 5 also shows values taken from the tabula- 
tions of Meisl and Shapiro: ° Makansi, Selke, and 
Bonilla;'® and the recent experimental work of 
Achener and Jouthas.”° Reference 25 recommended 
the relation 


AH,, = 2023 — 0.2129/ 


(10) 
1100 = ¢ = 1700°F 


At temperatures below about 1600°F, the quasi- 
chemical curve agrees more closely with that of 
Achener and Jouthas than does the virial curve. 
Because of this and the close agreement between the 
virial and quasichemical methods above 1600°F, 
Eq. 9 is recommended for computations of the heat of 
vaporization. 


ENTHALPY AND ENTROPY OF LIQUID AND VAPOR 


For saturated liquid the enthalpy and entropy are 
computed from Eqs. 11 and 12 


H!— (H3)° = 22.991(0.3893527 — 0.552955 
x 10-47? + 0.113726 x 10-'T3_— 20.023) (11) 


s! = 22.991(0.896497 log 7 — 1.10557 x 1047 
+ 0.170408 x 10~'7* — 1.792026) (12) 


For saturated vapor the corresponding properties 
are computed from the saturated liquid properties 
and the quasichemical heat of vaporization, Eq. 9. It 
can be shown however that, with the heat of vapor- 
ization defined by Eq. 9, the vapor properties are in 
fact independent of those of the liquid. The vapor 
composition is found by solving Eqs. 13 to 20; mean 
molecular weight is then obtained from Eq. 21. 





i (13) 
a NP 
N, 
ky= Nips (14) 
where 
Bz 
log Ry = A» + - (15) 
Bx 
log ky = —A, + =e (16) 
T 
and 
N,+N,+N,=1 (17) 
P3k,N} + Pk,N} + N,-1=0 (18) 


is solved numerically and values for N, and N, are 
computed. 


N, = kyPN? (19) 
Ny 1-N,+% (20) 

Thus 
M = M,(N, + 2Nz + 44) (21) 


which is used in Eq. 22 to compute the mass specific 
volume of the vapor 


V =RT/MP (22) 
where 


M ‘7 MN, si M2N2 = M,(1 + No) (23) 


If given the temperature of the saturated vapor, its 
composition, mean molecular weight, and specific 
volume are computed as shown above. Ewing”® gave 
the following equations for the enthalpy and entropy 
of the vapor in state ?: 


H*( T, P) = H°8(T) + X_2AHp + Xy AH (24) 
and 
ae cntn., RepSHe!. NGAHe) 
SF A = S°8 / é é4 — l P 
4 T, P) = S°8(T) a ee n 

_ RX: RX, 

+ 2 ln ky + =F Ink 

a 4M, seal 


_ RUN; In Ny + Nz INN2 + Ng In Ny) 


M (29) 





*Ewing discussed techniques for evaluating the constants 
A,, By, Ay, and B, from P—V-T data.'??6 From the data, 
Stone!? obtained A, = 4.3249, B, = 7204.2, A, = 10.6798, and 
B, = 16,325. 
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Equations 24 and 25 are used to compute the change 
in enthalpy and entropy along a constant-temperature 
path from saturated vapor at 7 and P, to an ideal 
monomeric gas and then to superheated vapor con- 
taining monomer, dimer, and tetramer. 

The enthalpy and entropy of the subcooled liquid is 
computed by starting with the corresponding proper- 
ties of the saturated liquid. The enthalpy change is 
given by 


dH, ==\1+ mas, dp (26) 
p, p, oT}, 
and the entropy change by 
; K /0P, 
dSp = ace dp (27) 
1 \? p 


As an approximation, (3p,/aT), is computed from 
Eq. 1. 


FORTRAN SUBROUTINE 


A FORTRAN subroutine was written, based on the 
quasichemical approach for computing sodium ther- 
modynamic properties in the range 500 to 2500°F. 
This subroutine was based on data from Refs. 12, 17, 
20, and 23. For a given temperature and pressure, 
the subroutine determines the state and computes the 
properties. The data of Ref. 17 were used to compute 
vapor pressure for temperatures =1600°F, and those 
of Ref. 12 for temperatures 2=1600°F. 

Table 4 is a listing of the FORTRAN subroutine. It 
can be used in any compatible computer program 
calling for the given properties, or as the core ofa 
table-making program. The values shown in Table 5 
are the thermodynamic and related properties of the 
saturated liquid and vapor as computed by the sub- 
routine of Table 4. 


SPECIFIC HEAT OF VAPOR 


The FORTRAN subroutine (Table 4) was used to 
compute the specific heat of the vapor at constant 
pressure along the saturation curve, taking 

(cf), = H*(T + 1,P,) — H&(7) (28) 
The values thus computed are sensitive to errors in 
the enthalpy formulation, and hence these results are 
listed in a separate table on a provisional basis (dis- 
cussed under “Transport Properties” which follows). 

The general thermodynamic approach presented 
here is tentative; it is based on the best experimental 
data available. However, if more high-quality experi- 
mental data were available, a different approach 
might be selected. That is, precise P—V—T mea- 
surements to temperatures well below 1700°F might 
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lead to the choice of true gas-base calculations, 
rather than the hybrid approach taken. On the other 
hand, extension of the precision saturated-liquid 
specific-heat measurements to temperatures well 
above 1600°F would probably strengthen the choice 
of true liquid-base calculations. 


Transport Properties 


The thermal conductivity and viscosity of both the 
saturated liquid and vapor phases have been re- 
viewed. For the range 210 to 2500°F, Table 6 sum- 
marizes the recommended values for the transport 
properties; Table 6 also lists liquid-phase Prandtl 
number (Pr), thermal diffusivity (@), specific heat at 
constant pressure (c,), electrical resistivity (p,), 
and surface tension (c). Table 7 contains provisional 
values of specific heat, Prandtl number, and thermal 
diffusivity for the saturated equilibrium vapor; there 
are relatively large discrepancies in the supporting 
data for Table 7. 


LIQUID THERMAL CONDUCTIVITY 


Nine different reviews and experiments (Fig. 6) 
were analyzed to determine the best relation for 
the thermal conductivity (Refs. 4, 5, 14, and 27 to 
32). In several instances thermal-diffusivity values 
were actually measured;”?:” they were converted to 
thermal-conductivity values using the specific-heat 
data of Ginnings.” We recommend Ginnings specific- 
heat relation, which, in engineering units, is 


c,[Btu/(1b)(°F)] = 0.389352 
~ 1.10599 x 10-77 + 3.41178 x 107°7? (29) 


where J is in °R, 
The analysis of these different reviews showed that 
the relation 


k[watts/(em)(°C)] = 0.916 — 0.49 x 107¥(°C) (30) 
of Evangelisti and Isacchini’’’® (Fig. 6) was in good 
agreement with other data from the melting point to 
1200°F. Values greater than 1200°F were computed 
from the Wiedeman-— Franz— Lorenz law using a fitted 
Lorenz value of 2.34 x 107° volts’/°K? and recom- 
mended electrical-resistivity values taken from 
Eq. 35. A single relation representing the thermal 
conductivity over the entire range from the melting 
point to 2500°F resulted from a smooth least-squares 
fit that was obtained with the use of values of the two 
proposed relations evaluated at 100°F temperature 
intervals from 200 to 2500°F. The result is given by 


k = 54.306 — 1.878 x 10~% + 2.0914 x 10~%? (31) 


where kand¢ are in Btu/(hr)(ft)(°F) and °F. 
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Table 6 TRANSPORT AND RELATED PROPERTIES OF SODIUM 



























: : ; Specific Thermal Electric 
Viscosity, Thermal conductivity, heat, Prandtl diffusivity, Surface resis- 
Temp., Ibm/(ft)(hr) Btu/(hr)(ft)(°F) liquid, number, liquid, tension, tivity, 
°F Liquid Vapor Liquid Vapor Btu/(lb,,)(°F) liquid sq ft/hr lb¢/ft Hohm-cm 
210.0 1.6657 0.03598 50.45 0.006285 0.010914 x 0.013489 9.67 
240.0 1.5374 0.03623 49.91 0.007347 0.010122 2.634 0.013375 10.29 
270.0 1.4272 0.03647 49.38 0.008391 0.009444 2.632 0.013261 10.91 
300.0° 1.3316 0.03672 48.86 0.009419 0.008858 2.630 0.013147 11.53 
330.0 1.2480 0.03696 48.33 0.010429 0.008347 2.627 0.013032 32.57 
360.0 1.1745 0.03721 47.81 0.011421 0.007899 2.623 0.012918 12.81 
390.0 1.1093 0.03745 47.29 0.012396 0.007505 2.619 0.012804 13.46 
420.0 1.0512 0.03770 46.78 0.013353 0.007155 2.615 0.012689 14.12 
450.0 0.9991 0.03794 46.27 0.014293 0.006843 2.610 0.012575 14.79 
480.0 0.9522 0.03819 45.77 0.015216 0.006564 2.604 0.012461 15.47 
510.0 0.9098 0.03843 45.27 0.016121 0.006314 2.5 0.012346 16.16 
540.0 0.8712 0.03868 44.77 0.017009 0.006088 2.6 0.012232 16.85 
570.0 0.8360 0.03893 44.28 0.017879 0.005884 0.012118 17.56 
600.0 0.8038 0.03917 43.79 0.018732 0.005699 0.012004 18.29 
630.0 0.7742 0.03942 43.30 0.019567 0.005531 0.011889 19.02 
660.0 0.7470 0.03966 42.82 0.020385 0.005377 2.562 0.011775 
690.0 0.7218 0.03991 42.34 0.021185 0.005238 2.553 0.011661 
720.0 0.6984 0.04015 41.86 0.021968 0.005110 54 0.011546 
750.0 0.6767 0.04040 41.39 0.022733 0.004993 0.011432 
780.0 0.6564 0.04064 40.93 0.023481 0.004886 a. 0.011318 
810.0 0.6375 0.04089 40.46 0.024212 0.004788 2.513 0.011203 
840.0 0.6198 0.04113 40.00 0.024925 0.004698 2.502 0.011089 
870.0 0.6033 0.04138 39.55 0.025621 0.004616 2.490 0.010975 
900.0 0.5877 0.04162 39.09 0.026299 0.004540 2.478 0.010861 
930.0 0.5730 0.04187 38.64 0.026960 0.004471 2.466 0.610746 
960.0 0.5592 0.04211 38.20 0.027603 0.004407 2.453 0.010632 28.00 
990.0 0.5461 0.04236 37.76 0.028229 0.004349 2.439 0.010518 28.91 
1020.0 0.5338 0.04260 37.32 0.028837 0.004296 2.425 0.010403 29.84 
1050.0 0.5221 0.04285 36.89 0.029428 0.004247 2.411 0.010289 30.79 
1080.0 0.5110 0.04310 36.46 0.030002 0.004203 2.396 0.010175 31.76 
1110.0 0.5004 0.04334 36.03 0.030558 0.2998 0.004163 2.381 0.010060 
1140.0 0.4904 0.04359 35.61 0.031096 0.2997 0.004127 366 0.009946 33. 
1170.0 0.4808 0.04383 35.19 0.031617 0.2997 0.004095 0.009832 34.78 
1200.0 0.4717 0.04408 34.78 0.032121 0.2997 0.004066 0.009718 35.83 
1230.0 0.4630 0.04432 34.37 0.032607 0.2998 0.004040 2.317 0.009603 36.90 
1260.0 0.4547 0.04457 33.96 0.033076 0.3000 0.004017 2.300 0.009489 37.99 
1290.0 0.4468 0.04481 33.56 0.033527 0.3002 0.003998 2.282 0.009375 39.10 
1320.0 0.4392 0.04506 33.16 0.033961 0.3005 0.003981 2.264 0.009260 40.24 
350.0 0.4319 0.04530 32.76 0.034377 0.3009 0.003967 2.246 0.009146 41.40 
1380.0 0.4249 0.04555 32.37 0.034776 0.3013 0.003955 2.228 0.009032 42.58 
1410.0 0.4181 0.04579 31.98 0.035158 0.3018 0.003946 2.209 0.008917 43.78 
1440.0 0.4117 0.04604 31.59 0.035522 0.3023 0.003939 2.190 0.008803 45.01 
1470.0 0.4055 0.04628 31.21 0.035868 0.3029 0.003935 2.172 0.008689 46.27 
1500.9 0.3995 0.04653 30.84 0.036197 0.3036 0.003933 2.152 0.008575 47.55 
1530.9 0.3937 0.04677 30.46 0.036509 0.3043 0.003933 2.132 0.008460 48.86 
1560.0 0.3882 0.04702 30.09 0.036803 0.3051 0.003935 2.112 0.008346 50.19 
1590.0 0.3828 0.04726 29.73 0.037080 0.3059 0.003940 2.092 0.008232 51.55 
1620.0 0.3777 0.04751 29.37 0.037339 0.3069 0.003946 2.071 0.008117 52.94 
1650.0 0.3727 0.04776 29.01 0.037581 0.3078 0.003955 2.051 0.008003 54.35 
1680.0 0.3678 0.04800 28.65 0.037805 0.3089 0.003965 2.030 0.007889 55.79 
1710.0 0.3632 0.04825 28.30 0.038012 0.3099 0.003977 2.009 0.007774 57.26 
1740.0 0.3587 0.04849 27.96 0.038202 0.3111 0.003991 1.988 0.007660 58.76 
1770.0 0.3543 0.04874 27.61 0.038374 0.3123 0.004007 1.967 0.007546 60.29 
1800.0 0.3501 0.04898 27.27 0.038528 0.3136 0.004025 1.946 0.007432 61.85 
1830.0 0.3459 0.04923 26.94 0.038665 0.3149 0.004044 1,925 0.007317 63.44 
1860.0 0.3420 0.04947 26.61 0.038785 0.3163 0.004066 1.904 0.007203 65.06 
1890.0 0.3381 0.04972 26.28 0.038887 0.3178 0.004089 1.882 0.007089 66.72 
1920.0 0.3343 0.04996 25.95 0.038972 0.3193 0.004114 1.861 0.006974 68.40 
1950.0 0.3307 0.05021 25.63 0.039039 0.3209 0.004140 1.839 0.006860 70.12 
1980.0 0.3272 0.05045 25.32 0.039089 0.3225 0.004168 1.818 0.006746 71.87 
2010.0 0.3237 0.05070 25.00 0.039121 0.3243 0.004198 1.796 0.006631 73.65 
2040.0 0.3204 0.05094 24.69 0.039136 ¢ 0.004230 1.775 0.006517 75.47 
2070.0 0.3171 0.05119 24.38 0.039134 0.004263 1.753 0.006403 77.32 
2100.0 0.3139 0.05143 24.09 0.039114 0.004298 1.732 0.006289 79.20 
2130.0 0.3109 0.05168 23.79 0.039076 0.004335 1.711 0.006174 81.12 
2160.0 0.3079 0.05193 2° 0.039021 0.004373 1.689 0.006060 83.08 
2190.0 0.3049 0.05217 2° 0.038949 0.004413 1.668 0.005946 85.07 
y 0.3021 0.05242 2 0.038859 0.004454 1.647 0.005831 87.10 
0.2993 0.05266 2 0.038752 0.004497 1.626 0.005717 89.17 
0.2966 0.05291 22. 0.038627 0.004542 1.605 0.005603 91.27 
0.2939 0.05315 22.08 0.038485 0.004589 1.585 0.005488 93.41 
0.2913 0.05340 21.81 0.038325 0.004637 1.564 0.005374 95.59 
0.2888 0.05364 21.54 0.038148 0.004686 1.544 0.005260 97.81 
0.2863 0.05389 21.28 0.037954 0.004738 1.523 0.005146 100.06 
0.2839 0.05413 21.02 0.037742 0.004791 1.503 0.005031 102.36 
0.2816 0.05438 20.76 0.037512 0.004845 1.483 0.004917 104.70 
0.2793 0.05462 20.51 0.037265 0.004901 1.463 0.004803 107.08 
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Table 7 PROVISIONAL VALUES FOR SPECIFIC HEAT, PRANDTL NUMBER, 
AND THERMAL DIFFUSIVITY OF SATL RATED SODIUM VAPOR 

































































Specific Thermal | Specific Thermal 
Temp., heat, Prandtl diffusivity, Temp., heat, Prandtl diffusivity, 
F Btu/(1b m) (CF) number sq ft/hr | F Btu/(1b»,)(°F) number sq ft/hr 
1 
500.0 0.349 0.846 0.387E 06 | 1550.0 0.628 0.804 0.483E 01 
530.0 0.366 0.846 0.196E 06 \| 1580.0 0.626 0.798 0.420E 01 
560.0 0.384 0.848 0.104E 06 | 1610.0 0.623 0.793 0.369E 01 
590.0 0.402 i 0.569E 05 | 1640.0 0.619 0.788 0.325E 01 
620.0 0.420 0.857 0.322E 05 | 1670.0 0.616 0.783 0.287E 01 
650.0 0.438 0.862 0.188E 05 1700.0 0.613 0.779 = 01 
680.0 0.456 0.868 0.112E 05 || 1730.0 0.610 0.775 > O01 
710.0 0.473 0.874 0.697E 04 | 1760.0 0.607 0.771 2E 01 
740.0 0.490 0.879 0.441E 04 \| 1790.0 0.604 0.768 E01 
770.0 0.506 0.884 0.286E 04 || 1820.0 0.601 0.766 E Ol 
800.0 0.888 0.189E 04 | 1850.0 0.599 0.763 E01 
830.0 0.892 0.128E 04 ] 1880.0 0.596 0.762 2E 01 
860.0 0.895 0.883E 03 | 1910.0 0.594 0.761 E01 
890.0 0.563 0.897 0.620E 03 1940.0 0.591 0.760 01 
920.0 0.574 0.898 0.442E 03 1970.0 0.589 0.760 E 00 
950.0 0.585 0.898 0.321E 2000.0 0.587 0.760 2E 00 
980.0 0.595 0.897 0.236E 2030.0 0.585 0.761 00 
1010.0 0.603 0.896 0.176E 2060.0 0.583 0.762 0.742E 00 
1040.0 0.611 0.894 0.133E 03 2090.0 0.581 0.763 0.678E 00 
1070.0 0.618 0.891 0.102E 03 2120.0 0.580 0.765 0.621E 00 
1100.0 0.623 0.888 0.794E 02 2150.0 0.578 0.768 0.569E 00 
1130.0 0.628 0.884 0.622E 02 2180.0 0.577 0.771 0.523E 00 
1160.0 0.632 0.879 0.493E 02 2210.0 0.575 0.774 0.481E 00 
1190.0 0.635 0.874 0.394E 02 2240.0 0.574 0.443E 00 
1220.0 0.637 0.869 0.318E 02 2270.0 0.572 0 2 0.408E 00 
1250.0 0.639 0.863 0.259E 02 25 0 0.571 0.786 0 > 00 
1280.0 0.640 0.858 0.212E 02 25 0 0.569 0.791 0 > 00 
1310.0 0.640 0.852 0.175E 02 2 0 0.568 0.796 0.322E 00 
1340.0 0.640 0.845 0.146E 02 2 0 0.566 0.801 0.299E 00 
1370.0 0.639 0.839 0.122E 02 2 .0 0.564 0.807 0.277E 00 
1400.0 0.638 0.833 0.103E02 || 2450.0 0.563 0.813  0.257E 00 
1430.0 0.637 0.827 0.880E 01 2480.0 0.561 0.819 0.239E 00 
1460.0 0.635 0.821 0.752E 01 
1490.0 0.633 0.815 0.645E 01 
1520.0 0.631 0.809 0.557E 01 \| 
60 | 
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Fig. 6 Comparison of analyses for the 


rl 5 = ssh ayes, 
thermal conductivity of liquid sodium. 
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‘ Dashed lines are extrapolated or calculated 
\y values. Recommended curve of Et angelisti 
and Isacchini?”’ includes a high-tempera- 
10 | | | | ture portion computed from the Wiedeman 
200 500 1000 1500 : 2000 2500 Franz —Lorenz law. 
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VAPOR THERMAL CONDUCTIVITY 


To our knowledge the only experimental thermal- 
conductivity values published for sodium vapor are 
those of Stefanov.*® The measurements were made 
with a dilatometer described by Timrot and Totskii.* 


A number of experimental and theoretical studies 


of the conductivity values given in Ref. 33 to obtain 
values for Tables 6 and 7. The result of this fit is 


k = 0.1639 x 107 + 0.3977 
x 10~4¢ — 0.9697 x 10-84? = (32) 


where k and ¢ are in Btu/(hr)(ft)(°F) and °F. 


have been conducted (Refs. 7, 8, 35, and 36) or are 
under way *’-“ to determine vapor conductivity. Fig- 


ure 7 shows the results reported. There is little LIQUID VISCOSITY 


The viscosity of liquid sodium has been measured 

| l from the melting point to 1920°F by the investigators 

a shown in Table 8. The methods used to measure the 

viscosity, the equation, and references are also listed 
in Table 8. 

Several viscosity compilations have appeared in 
literature, and some are listed in Table V.3of Ref. 1. 
The correlations of these data show reasonably good 
agreement between 212 and 1650°F. Above 1830°F 
the deviations become larger. 

To obtain a relation for the liquid viscosity over 
the range from melting point to 2500°F required that 
a smooth curve be drawn through a plot of all data 
points. The best fit to the curve is given by 


log n = 1.0203 + a — 0.4925 log T (33) 





] 


Lee & Bonilla‘? 
Stefanov 33 (exp.) 














500 1000 


1500 2000 2500 


Fig. 7 Curves of thermal conductivily of sodium vapor 
show little agreement. Data of Stefanov® are arbitrarily 
recommended until the work of Lee*” becomes available. 


where 7 and 7 are given in lb,,/(hr)(ft) and °R. A plot 
of Eq. 33 is shown in Fig. 8; calculated values are 
agreement among most of the sources; in fact, the given in Table 6. 
predictions of Lee and Bonilla“’ and the experimental 
results of Stefanov’® disagree as much as 25%. How- 
ever, together their results are the best information 
available. Arbitrarily we recommend the data of 
Stefanov, at least until the Lee’ results become 
available. A second-degree polynomial fit was made 


VAPOR VISCOSITY 
Stefanov reported smoothed experimental results 
for the temperature range 1263 to 1598°F and pres- 
sures 0.055 to 0.398 kg/cm’, A linear least-squares 








Table 8 LIQUID-SODIUM VISCOSITY INVESTIGATIONS 
Temp. 
Investigators Method range, °F Equation Refs. 
Chiong™ Oscillating sphere 208 to 672 *nv's = 1.183 x 10~3e M6 -5/eT Hagen"® 
Godfrey” Oscillating sphere 620 to 1612 T 
Ewing®! Modified Ostwald 218 to 379 *nu's = 1.089 x 10~3¢ 756 -8/eT Ewing'® 
glass viscometer 

Ewing*! Nickel viscometer 289 to 1267 *nu’ = 1.142 x 10~3e 739-8/uT Ewing!® 
(Ostwald type) 

Novikov®:83 Oscillating cylindrical 212 to 1470 t Shvidkovsky ®4 
capsule Stone!2 

Shpilrain*4 Oscillating cylinder 320 to 1470 *log n = —1.6814 — 0.4296 log T + 234.65/T —Goltsova!’ 


Achener® Oscillating cylinder 208 to 1920 flog n = 3.25538 — 0.96843 log T Shvidkovsky *4 


_ 10.4126 
ra 
or 
tnv's = 1.54222 e506 -151/vT 


All the above 
references 


290 65 
Golden and Tokar! Compilation M.p. to 2500 Slog 7 = 0.5108 + a — 0.4925 log T 





*7 is the viscosity in poises, v is the specific volume in cm?/g, and T is in degrees Kelvin. 
TData presented in Table V.2 in Ref. 1. 

tn is viscosity in millipoises. 

§n is in centipoises. 


REACTOR AND FUEL-PROCESSING TECHNOLOGY, Vol. 11, No. 1, Winter 1967-1968 





FLUID AND THERMAL TECHNOLOGY 39 


































T 7 a a. ae T 
Stenfanov’’(calc.; monomer vapor) 
Weatherford®(calc.; sat. equil. vap.) 
Lory - Stefanov33(exp.: sat. equil. vap.) 
— . (6Fr- 4 
: < 
| oo : _ & 
= 10 Golden & Tokar * sa 
2 P= 
2 D 
S ° 
2 3 
> =~ 
eo = 
Weatherford®”(calc.; sat. equil. vap 
| | | | | me 3 | 
500 1000 1500 2000 2500 500 1000 2000 2500 
Temperature, °F Temperature. °F 
Fig. 8 The viscosity of liquid sodium. Fig. 9 Curves of viscosity of sodium vapor show Little 


agreement. Data of Stefanov® are presently recommended. 
fit was made of the already smoothed values to 


; recommend that the Stefanov data (Eq.34) be used 
result in 


until additional experimental information becomes 

= ee available. Calculated values for the vapor viscosity 

Ib,,/(hr)(ft)] = 0.03427 + 8.176 x 10~%4(° 
{bm / (hr) (ft) . = ae) oH are given in Table 6. 

The curve calculated from Eq. 34 is plotted in Fig. 9 

along with other curves derived by experimental and 


sie he nets : Electrical Resistivity 
theoretical investigations. Table 9 lists the investi- 


gators, the method used, and the references. 
Because there is moderate agreement among most 
of the recent theoretical and experimental efforts, we 


The electrical-resistivity values from Tepper” 
cover the widest experimental temperature range: 
they also represent an average of the divergent ex- 


trapolated values from other investigators (Table 10). 
A least-squares fit to the Tepper data is given by 





Ref. 2.1737 + 7.6254 x 107% 








) 
Investigators Method sources Pe ‘ a 
ea + 5.8207 x 10° ¢° + 1.1262 x 10 “#" = (35) 
Lee and Bonilla*® Experimental 87 
y (preliminary) ; 
Dunning’ a, Calculated 88 where p, is in pohm-in. and / is in F. 
Weatherford® Calculated Not clear T 1 ical istivity 1 fr fi f th 
Weatherford® Calculated 39 he electrical-resistivity values from five o e 
Davies” Calculated None referenced investigators are in very good agreement 
anov® Experimental None * 
— fact ee with the values from Tepper for the range 300 to 
Table 10 ELECTRICAL-RESISTIVITY INVESTIGATIONS OF LIQUID SODIUM 





Temp. 
Investigators Method range, °F Equation 





Evangelisti and Isacchini2’+%° Parallel resistor 212 to 1832 *p, = 6.56 + 2.868 x 10-2 

+ 1.868 x 10742 + 9.078 x 10-%3 
Tepper!4 Parallel resistor 272 to 1989 tp, = 2.1729 + 7.6248 x 10 3t : 

+ 5.8313 x 1077t2 + 1.1260 x 107%¢ 
Kapelner and Bratton*! Parallel resistor 208.4 to 1567 *p, = 5.904 + 3.695 x 10°-*t — 9.359 

x 10° "¢? + 2.291 x 107%? 
Freedman and Robertson*! Parallel resistor 212 to 1022 *4y e ~ 6.3715 + 2.8725 « 10 2t + 1.8396 

x 107542 + 5.1593 x 107%3 
Rahiser*? Parallel resistor 392 to 1472 *ip, = 7.6986 + 2.3705 x 10 2t 

+ 3.0095 x 107-%¢2 
Semyachkin and Solovev*? Parallel resistor 212 to 1832 *tp, 6.4525 + 3.3775 x 10° + 9.127 

x 10782 + 1.12515 x 107843 





*p, is inohm-cm and ¢ is in °C. 
tp, is in wohm-in. and ¢ is in °F. 
+Curves were fitted to investigators’ data. 
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Fig. 10 Recommended electrical-resislivily values for 


liquid sodium are based on data by Tepper."4 


1600°F, Also, as stated before, Eq. 35 represents an 
average of the extrapolated values for temperatures 
outside the area of agreement. It is recommended 
that the Tepper relation be used with temperatures 
from the melting point to 2500°F. Figure 10 plots 
the recommended Eq. 35; Table 6 lists values of 
electrical resistivity for various temperatures. 


Surface Tension 


The surface-tension data are shown in Fig. 11. 
Except for the Poindexter and Kernaghan‘! curve, 
there is moderately good agreement among all the 
correlations both with respect to the temperature 
coefficient and the actual surface-tension values 
obtained near the melting point (Table 11). As a 
result, choosing a recommended fit was difficult. So 
the data were pooled and the values averaged for 
surface tension at the melting point and the tempera- 
ture coefficients (Poindexter and Kernaghan results 
were excluded). The resulting average melting-point 






















| | 
Poindexter & Kernaghan 4! 
| a 
Addison 
. Jordan & Lane9? 
5 150  Addison% - 
: Taylor99-101 Solovev & Makarova 
= 
8 Tepper!02-104 
F 100 _ 
Golden & Tokar! —— 
60 | | 
200 500 1000 1500 2000 2500 
Temperature, °F 
Fig. 11 Surface-tension data for sodium are in good 


agreement; the recommended curve (Golden and Tokar") is 
an average of various referenced data, 
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surface-tension value and temperature coefficient 
are 196.9 dynes/cm and —0.10 dyne/(cm)(°C). The 
recommended relation of surface tension and tem- 
perature is given by 


o = 206.7 — 0.10/ (36) 


where o and / are in dynes/em and °C. Equation 36 
can be rewritten as 


o = 0.01429 — 3.81 x 107° (37) 


where o and / are in units of lb,/ft and °F. Surface- 
tension values are given in Table 6. 








Table 11 SURFACE-TENSION DATA OBTAINED 
FROM CORRELATIONS 
Surface Temp. 
Temp. tension ,* coefficient, 
Investigators range, °C  dynes/em_ dynes/(cm)(°C) 
Addison*4 128 to 180 195 —0.10 
Addison®*:*6 110 to 220 202 —0.10 
Jordan and Lane*’ 130 to 245 200.2 —0.11 + 0.01 
Poindexter and 
Kernaghan"! 103 to 243 206.5 —0.046 

Solovev and Makarova*’ 194 to 932 200.7 —0.0894 
Taylor®?-1% 110 to 450 190.8 —0.10 
Tepper102- 104 141 to 992 192.5 —0.0986 





*At melting point, 98°C. 


Thermal-Expansion Coefficient 


Table 12 gives values for the thermal-expansion 
coefficient, together with related data on compress- 
ibility and sonic velocity. In the absence of direct 
dilatometric measurements of thermal expansion at 
high temperatures, these values of the thermal- 
expansion coefficient, @,, have been estimated from 
Eq. 1 as follows: 


“4 .) or ee a) 
a, = 767 oe =H(35 , -” 


where the subscript s indicates saturation conditions. 
Figure 12 shows the thermal-expansion curve com- 
puted from this approximation. 

Equation 38 is valid only for liquid sodium in 
equilibrium with its saturated vapor; a small cor- 
rection can be made if the pressure is significantly 
greater than the vapor pressure. For this condition 
the isothermal compressibility 


__1/av) _ 1/ep 
B, = = al (39) 


is used in the equation 


gi Farien. a+ B,(57) sii 


p? 
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Table 12 THERMAL-EXPANSION, COMPRESSIBILITY, AND THERMAL-PRESSURE 
COEFFICIENTS FOR LIQUID SODIUM 





























Thermal- isothermal Thermal- Specific 
expansion compress- pressure Sonic Adiabatic heat, 
Temp., coefficient, ibility, coefficient, velocity, compressibility, const. vol., 
°F by atm“! atm/°F ft/sec atm"! Btu/(1b »)(°F) 
210 1.3930-004* 1.7663-005* 7.89 8285 1.7150-005 0.3573 
240 1.4014-004 1.7888-005 7.83 8257 1.7342-005 0.3559 
270 1.4098-004 1.8117-005 7.78 8228 1.7537-005 0.3545 
300 1.4183-004 1.8350-005 7.73 8199 1.7734-005 0.3530 
330 1.4268-004 1.8587-005 7.68 8171 1.7935-005 0.3516 
360 1.4353-004 1.8829-005 7.62 8142 1.8140-005 0.3503 
390 1.4438-004 1.9075-005 7.57 8113 1.8347-005 0.3489 
420 1.4524-004 1.9325-005 7.52 8085 1.8558-005 0.3475 
450 1.4610-004 1.9580-005 7.46 8056 1.8772-005 0.3461 
480 1.4696-004 1.9840-005 7.41 8027 1.8989-005 0.3447 
510 1.4783-004 2.0105-005 7.35 7999 1.9210-005 0.3433 
540 1.4870-004 2.0374-005 7.30 7970 1.9435-005 0.3420 
570 1.4957-004 2.0649-005 7.24 7941 1.9663-005 0.3406 
600 1.5044-004 2.0928-005 7.19 7913 1.9895-005 
630 1.5132-004 2.1213-005 7.13 7884 2.0131-005 
660 1.5219-004 2.1503-005 7.08 7855 2.0371-005 
690 1.5308-004 2.1799-005 7.02 7827 2.0614-005 
720 1.5396-004 2.2100-005 6.97 7798 2.0862-005 
750 1.5485-004 2.2406-005 6.91 7769 2.1114-005 
780 1.5574-004 2.2719-005 6.86 7741 2.1370-005 
810 1.5663-004 2.3037-005 6.80 7712 2.1630-005 
840 1.5753-004 362-005 6.74 7683 2.1894-005 
870 1.5843-004 2.3693-005 6.69 7655 2.2163-005 
900 1.5933-004 2.4030-005 6.63 7626 2.2437-005 
930 1.6024-004 2.4373-005 6.57 7598 2.2715-005 
960 1.6115-004 2.4723-005 6.52 7569 2.2998-005 0.32¢ 
990 1.6206-004 2.5080-005 6.46 7540 2.3286-005 0.3219 
1020 1.6298-004 2.5443-005 6.41 7512 2.3578-005 0.3206 
1050 1.6389-004 2.5814-005 6.35 7483 2.3876-005 0.3193 
1080 1.6482-004 2.6191-005 6.29 7454 2.4179-005 0.3180 
1110 1.6574-004 2.6576-005 6.24 7426 2.4487-005 0.3168 
1140 1.6667-004 2.6969-005 6.18 7397 2.4800-005 0.3155 
1170 1.6760-004 2.7369-005 6.12 7368 2.5119-005 0.3142 
1200 1.6854-004 2.7777-005 6.07 7340 2.5444-005 0.3130 
1230 1.6947-004 2.8193-005 6.01 7311 2.5774-005 0.3117 
1260 1.7042-004 2.8618-005 5.95 7282 2.6110-005 0.3104 
1290 1.7136-004 2.9050-005 5.90 7254 2.6452-005 0.3092 
1326 1.7231-004 2.9491-005 5.84 7225 2.6800-005 0.3080 
1350 1.7326-004 2.9941-005 5.79 7196 2.7154-005 0.3067 
1380 1.7422-004 3.0400-005 5.73 7168 2.7515-005 0.3055 
1410 1.7517-004 3.0868-005 5.68 7139 2.7882-005 0.3043 
1440 1.7614-004 3.1345-005 5.62 7110 2.8255-005 0.3031 
1470 1.7710-004 3 5.56 7082 2.8636-005 0.3019 
1930 1.9237-004 4.0662-005 4.73 6642 3.5435-005 
1960 1.9340-004 4.1340-005 4.68 6614 3.5950-005 
1990 1.9443-004 4.2032-005 4.63 6585 3.6475-005 
2020 1.9546-004 4.2738-005 4.57 6556 3.7010-005 
2050 1.9650-004 4.3459-005 4.52 6528 3.7556-005 
2080 1.9754-004 4.4195-005 4.47 6499 3.8113-005 
2110 1.9859-004 4.4947-005 4.42 6471 3.8680-005 
2140 1.9964-004 4.5715-005 4.37 6442 3.9259-005 
2170 2.0069-004 4.6499-005 4.32 6413 3.9849-005 
2200 2.0175-005 4.7300-005 4.27 6385 4.0451-005 
2230 2.0281-004 4.8118-005 4.21 6356 4.1065-005 
2260 2.0388-004 4.8954-005 4.16 6327 4.1691-005 
2290 2.0495-004 4.9807-005 4.11 6299 4.2330-005 
2320 2.0602-004 5.0679-005 4.07 6270 4.2982-005 
2350 2.0710-004 5.1569-005 4.02 6241 4.3647-005 
1520 1.7872-004 3.2665-005 5.47 7034 2.9285-005 
1550 1.7970-004 3.3178-005 5.42 7005 2.9684-005 
1580 1.8067-004 3.3702-005 5.36 6977 3.0091-005 
1610 1.8166-004 3.4236-005 5.31 6948 3.0505-005 
1640 1.8264-004 3.4781-005 5.25 6919 3.0927-005 
1670 1.8363-004 3.5337-005 5.20 6891 3.1356-005 
1700 1.8463-004 3.5905-005 5.14 6862 3.1794-005 
1730 1.8562-004 3.6484-005 5.09 6833 3.2240-005 
1760 1.8663-004 3.7075-005 5.03 6805 3.2694-005 
1790 1.8763-004 3.7678-005 4.98 6776 3.3156-005 
1820 1.8864-004 3.8294-005 4.93 6747 3.3628-005 
1850 1.8965-004 3.8922-005 4.87 6719 3.4108-005 
1880 1.9067-004 3.9564-005 4.82 6690 3.4598-005 
1910 1.9169-004 4.0218-005 4.77 6662 3.5097-005 
2370 2.0782-004 5.2174-005 3.98 6222 4.4098-005 
2400 2.0890-004 5.3097-005 3.93 6194 4.4786-005 
2430 2.0999-004 5.4040-005 3.89 6165 4.5488-005 0.2676 
2460 2.1108-004 5.5003-005 3.84 6136 4.6205-005 0.2666 
2490 2.1218-004 5.5988-005 3.79 6108 4.6937-005 0.2657 
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Adiabatic compressibility is determined from the - 
o velocity of sound, C, by 
+ Bs = 1/pc? (43) 
t=] 
8 and p from Eq. 1. 
2 SONIC VELOCITY 
“ Sonic velocities in sodium have been measured by 
= a number of investigators, as indicated in Table 13. 
= - | | | | Pochapsky used triply distilled sodium and an 





200 500 1000 1500 2000 2500 acoustic interferometer to measure the temperature 


Temperature, °F 





Fig. 12. Estimated thermal-expansion coefficient for 


sodium. 


The derivative (dp/eT). can be obtained from the 
slope of the vapor-pressure equation using the equa- 
tions of Ditchburn and Gilmour?’ and of Stone" in the 
appropriate temperature range. 


Compressibility 


Although there have been many studies of the 
compressibility of solid sodium, no direct measure- 
ments of isothermal compressibility, /,;, have been 
made over the temperature range of interest, al- 1.6 | | | L | 
though Bridgeman” measured the change in com- 200 500 1000 1500 2000 2500 
pressibility at the melting point. The isothermal 
compressibility is related to the easily and accu- 
rately measured adiabatic compressibility, 6., by the 


Isothermal compressibility, B, 1079 /atm 














Temperature, °F 


Fig. 13 Calculated isothermal compressibility of sodium. 





relations 
Cy 
By =— Bs (41) ox . ; 43 
C, coefficient with excellent accuracy and consistency. 
and His temperature coefficient has been substantiated by 
| me. . Te others,“ including the most recent work of Trelin.“® 
Py * Ps * ae re © pcp (42) Therefore a linear equation, based on Pochapsky’s 


constants, for sonic velocity and temperature was 


where the subscript S refers to an adiabatic process. 
Values of the isothermal compressibility are plotted 
in Fig. 13. 


used with the density relation (Eq. 1) to calculate the 
adiabatic compressibility as a function of tempera- 
ture. 


Table 135 EXPERIMENTAL DETERMINATIONS OF SONIC VELOCITY IN LIQUID SODIUM 





Temp. 


Sonic velocity at coefficient, Temp. 





Investigators melting point, m/sec m/(sec)(°C) range, °C Comments 
Kleppa!® 2395 + 25 at 98°C =0.3 98 to 235 1% calcium and other impurities including 
gas 
Pochapsky*? 2526 + 5 at 97.6 + 0.1°C —0.524 + 0.003 97.6 to 272 Acoustic interferometer; good purity 
Abowitz and Gordon*4 2526 + 5 —0.524 + 0.003 100 to 140 
Trelin'™ 2595 + 47 —0.577 100 to 700 Pulsed interferometer 
Jarzynski and Litovitz!" 2526 —0.50 104 to 154 
Letcher and Beyer!" 2525 + 25 —0.50 97.8 to 150 
Ilgunas and Barsauskas!% 2533 —0.66 100 to 180 
Ying and Scott* 2526 at 97.8°C — 0.523 97.8 to 315 
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THERMAL-PRESSURE COEFFICIENT 


The heat capacity at constant volume, c,, and the 
thermal-pressure coefficient, y,,, where c, can be 
obtained from Eq. 41 and 


_ [ap _ Ap 
Vy =(5 7), 3, (44) 





were calculated. The thermal-pressure coefficient, 
which is plotted in Fig. 14, can be used to calculate 
the pressure change in confined sodium when heat is 
added from the relation 





>= (2) ar (2%) av 
dF -(5 7), dT + av ig (45) 
This gives 
i 1 AV 
AP = Ve AF = a (46) 


on substitution and integration. 





” 7 


Thermal — Pressure Coefficient, ¥,, atm/°F 














—— 
. _ 
= l l | : 
200 500 1000 1500 2000 2500 


Temperature, °F 


Fig.14 Calculated thermal-pressure coefficient of sodium. 


Summary 


The best values for the thermophysical properties 
of sodium, as presented here, are based on critical 
evaluations of the literature through 1966. Enthalpies 
and entropies (Figs. 15 and 16) were computed via a 
FORTRAN subroutine, as were many related factors. 
The recommended values have been presented pri- 
marily in tabular form to facilitate engineering use. 


Nomenclature 


C  =Sonic velocity; ft/sec 
Cp = Specific heat at constant pressure; Btu/(lb)(° F) 


c, = Specific heat at constant volume; Btu/(lb)(° F) 
AD}= Heat of dimerization 

E = Internal energy; Btu/lb, Btu/(lb)(mole) 

F = Free energy; Btu/lb, Btu/(lb)(mole) 

H =€Enthalpy; Btu/lb, Btu/(1b)(mole) 

K = R’'/R= 2.72131, Btu/(atm)(cu ft) 


k =Thermalconductivity; Btu/(hr)(ft)(° F), cal/ 
(sec)(em)(°C) 

M = Molecular weight (M; = 22.991) 

N = Mole fraction 

P = Pressure; atm abs. 


Pr = Prandtl number 1c)/k 

R =Gas constant; 0.730229 (atm)(cu ft)/(lb)(mole) 
(°R) 

R’ = Gas constant; 1.98718 Btu/(lb)(mole)(°R) 

S = Entropy; Btu/(lb)(°R), Btu/(lb)(mole)(°R) 


= Absolute temperature; °R(=°F + 459.7), 

°K(=°C + 273.15) 

t = Temperature; °F, °C 

V = Specific volume; cu ft/lb 

v  =Specific volume; cu ft/(lb)(mole)(monomer), 
em?/g 

xX = Weight fraction 

a = Thermal diffusivity; k/p,c,; sq ft/hr 

a, = Thermal-expansion coefficient; — 

@, = Thermal-expansion coefficient for saturation 
conditions; °F~' 

8, = Adiabatic compressibility; atm™' 

8B, = Isothermal compressibility; atm~! 

, = Thermal-pressure coefficient; atm/°F 

n = Dynamic viscosity; 1b,,/(hr) (ft) 

p = Vapor density; lb/cu ft, (1b)(mole)(monomer)/ 
cu ft 

p, = Liquid density; lb/cu ft, (1b)(mole)(monomer)/ 
cu ft 

p, = Electrical resistivity; .ohm-cm 

o =Surface tension; 1b;/ft 


SUBSCRIPTS 


D = Dimerization 

f = Formation 

i = State 

1 = Liquid 

m = Mass 

Pb = Pressure 

Ss = Adiabatic process; saturation 
T = Tetramerization 
l = Vaporization 

0 = Absolute zero 

1 = Monomer 

2 = Dimer 

4 = Tetramer 


SUPERSCRIPTS 
3 Standard state 
l = Liquid phase 
g Gaseous phase 
s  =Solid phase 


ll 


ll 
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AQUEOUS PROCESSING 


Research and Development 
on Aqueous Processing 


By C. E. Stevenson and D. M. Paige 


Caustic decladding followed by HNO, dissolution is 
favored as the feed-preparation method for process- 
ing D2O0-moderated organic-cooled reactor (HWOCR) 
fuels consisting of oxides or carbides clad in sin- 
tered aluminum product (SAP). A preliminary study 
indicates the feasibility of processing stainless- 
steel-clad oxide fuels from a fast breeder reactor 
in an existing commercial TBP-process_ solvent 
extraction plant. The chemistry and some process 
methods for transuranium elements are reviewed. 


Dissolution of Fuels Clad 
with Sintered Aluminum Product 


A recent report! from the Oak Ridge National 
Laboratory (ORNL) has reviewed potential methods 
of processing fueis proposed for use in a heavy- 
water-moderated organic-cooled reactor. Of the fuel 
materials being considered, carbides or oxides of 
uranium (or uranium and thorium) clad in SAP 
appear to be the leading candidates. Work at ORNL 
on the preparation of these fuels for processing was 
recently discussed in Power Reactor Technol. and 
Reactor Fuel Process.* The current report! indicates 
that caustic decladding of these fuels, followed by 
HNO, leaching, appears to be the most promising 
head-end process. Nitric acid dissolution of the 
cladding was very slow; only 18% of a sample of 
Alcoa XAP 001 material (2.4% O, 1.0% Si + Fe) was 
dissolved in 4 hr in boiling 8.7 HNO |average rate, 
0.17 mg/(min)(cm’) |, and catalysis by Hg’* or F7 did 
not markedly accelerate the rates of dissolution. It 
was of interest to note, however, that, when SAP was 
first heated to a temperature above its melting point, 
relatively rapid dissolution of the material’ was 
effected in 8M HNO3. Rapid dissolution of SAP also 
oecurred in boiling 3M to 12M HCl. A residue of 
Al,O; remained after treatment with acids or caustic. 


In NaOH or NaOH—NaNO; solutions (the presence 
of NaNO, suppresses evolution of hydrogen and yields 
NH; in its place), the aluminum present in SAP 
dissolves rapidly. Following removal of the caustic 
solution, the residual UO, or UC may be dissolved in 
HNO, (unirradiated UC appears to be hydrolyzed by 
caustic, whereas irradiated UC probably is not). In 
the HNO; dissolution the carbon present as carbide 
is incompletely oxidized to CO,; the balance is con- 
verted to organic acids, which may be extracted from 
strong HNO; solutions before uranium or uranium 
thorium extraction by TBP. A feed-preparation flow 
sheet is shown in Fig. 1. 
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Fig. 1 Caustic decladding of SAP-clad UC fuel.4 


Decladding of SAP-clad or zirconium-clad fuels 
by fluid-bed hydrochlorination, shearing and leaching, 
or NH,F treatment of zirconium cladding is also 
reviewed in application to HWOCR fuels, but these 
processes appear to be less promising. 
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Aqueous Processing 
of Fast Breeder Reactor Fuels 


Nicholson, Oak Ridge National Laboratory, has 
made a preliminary evaluation of the capability of 
the Purex processing facility of Nuclear Fuel Ser- 
vices, Inc., at West Valley, N. Y., for the processing 
of a fast breeder power-reactor fuel.’ For this pur- 
pose he has assumed a plant feed consisting of 
stainless-steel-clad UO,—PuO, fuel combined with 
the axial blanket. This plant feed, with atotal fission- 
product content equivalent to an average burnup of 
38,500 Mwd/metric ton and an average plutonium 
content equal to 10% of the heavy metals, corresponds 
to the fuel specification of a reactor design made by 
the General Electric Company.‘ A suggested flow 
sheet for the first-cycle extraction (codecontamina- 
tion and partition) is shown in Fig. 2. 


A daily throughput of 0.52 ton of fuel cooled for 
150 days was selected as being compatible with the 
existing plant operating license and within the capac- 
ity of the solvent extraction columns for the specified 
flow sheet. The shear—leach head-end facility was 
considered adequate for the previous fuel charge 
(with its increased heat dissipation and fissile con- 
tent) without major modifications. However, a new 
tail-end facility for plutonium purification and pack- 
aging was specified to achieve a plutonium through- 
put as high as required for the selected fuel-process- 
ing rate. This tail-end facility would consist of a 
shielded, continuous, anion-exchange system. Criti- 



































employed both concentration limits and geometrical 
restrictions for various parts of the plant. 

Calculations indicating a solvent irradiation dose 
(per pass) of 0.017 watt-hr/liter for fuel cooled 150 
days and 0.043 watt-hr/liter for fuel cooled 30 days 
were taken as a basis for the conclusion that solvent 
damage would be tolerable with the existing system 
at either cooling time, 

The present plant license for the release of Skr, 
‘317 and °H may be limiting for the processing of 
such quantities of fast breeder reactor fuel, partic- 
ularly if short decay cooling periods are involved, 
but Nicholson indicates that better performance of 
the off-gas treating equipment can be expected than 
is assumed in the license application and that greater 
allowances may be made for dilution. It therefore 
appears that an upward revision of the plant capacity 
might result if monitoring verifies the safety-analysis 
calculations. 


Processing of the Actinide Elements 


An excellent review article® by Los Alamos work- 
ers on the chemistry, preparation, and processing of 
the actinide elements has recently appeared. Produc- 
tion of the higher transuranium elements in the High 
Flux Isotope Reactor (HFIR) and their recovery in 
the Transuranium Facility (TRU) at Oak Ridge are 
discussed. The article also predicts that by the mid 
70’s we should see the production of about 1 g of 
*52Cf per year and of significant quantities of other 
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Fig. 2 
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First-cycle flow sheet for NFS processing of fast breeder reactor fuel.’ 
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Watanabe,® Japan Atomic Energy Research Insti- 
tute, has summarized some recent work of that 
Institute on the solvent extraction of heavy elements 
ranging from thorium to curium by tributyl phos- 
phine oxide (TBPO) and trioctyl phosphine oxide 
(TOPO). Distribution-ratio data for both the chloride 
and nitrate systems, summarized by Watanabe, are 
shown in Fig. 3. This figure includes new data ob- 
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Fig. 3 Extraction of actinide elements from HCl and HNO 


by TBP, TBPO, and TOPO.® 


tained at the Institute. The data indicate that TBPO 
may be used for separation of thorium or Np(V) 
from protactinium, U(VI), Np(IV), or Np(VI) in an HCl 
system. In the TOPO—HCI1 system, trivalent actinides 
may be separated from actinides having higher va- 
lences, but the separation of tetra-, penta-, and 
hexavalent species from each other is difficult. 
Uranium-VI may be separated from thorium and 
Np(V) by extracting from 1N to 2N HCl. In the TOPO-— 
HNO; system, americium and curium are separated 
from tetra- and hexavalent actinides and protac- 
tinium, and selective separations of the latter from 
elements other than actinides are also possible. 
These data were applied to effect various separa- 
tions. A solvent composed of 1% TOPO in toluene 
was used to extract 7°*U and 7**Pa from a solution of 
irradiated thorium dissolved in 6N HCl, and **’Pa 
was scrubbed from the extract with 6N HCl saturated 
with NaF. Americium was separated from rare earths 
(promethium, samarium, europium, and lanthanum) 
by reversed phase chromatography using columns 


containing TOPO or di-(2-ethylhexyl) phosphoric acid 
(HDEHP) absorbed on Kieselguhr and HNO; elution. 
Although either solvent could be used for separation 
of lanthanum, HDEHP was found to be better than 
TOPO for separating americium from promethium 
and europium. 


Plant and Equipment 
Design and Operation 


HYDROCLONES 


The work on stacked-clone liquid— liquid contactors 
has been reviewed previously.’-!! However, two re- 
cent reports are worthy of review since they cover 
work on a full-scale plant prototype model as well 
as data on 1.4-scale and 0.5-scale models. 

Whatley and Groenier!* of ORNL discuss the design 
features and present data on the full-scale (Mark-X 
size) 18-stage model. This prototype hydroclone is 
an integral-pump stacked-clone stainless-steel model 
(Figs. 4 to 7). Figure 8 is a sketch of the assembled 
18-unit machine as tested. They also present data on 
the 0.5-scale plastic model. 

Efficiencies for the prototype were compared with 
those for the 0.5-scale plastic model by using uranyl 
nitrate as the distribution species. The aqueous 
phase (continuous) was 1.7 NaNO;—0.08.\7 HNOs, and 
the organic phase was 18% TBP in Amsco. The 
authors report 52 to 82% efficiency for the prototype 
and 92 to 95% efficiency for the 0.5-scale model. 

A general review of hydroclones is also presented 
in the report. '” 

The effect of geometric scale-up on the perfor- 
mance of stacked-clone contactors is discussed ina 
report by Finsterwalder, '* working at ORNL as a 
guest scientist from West Germany. He presents 
data on a 1.4-scale model as well as on a 0.5-scale 
model and compares these with data on the full-scale 
Mark-X model. Table 1 presents a summary of the 
data for the three units. 

The hydraulic K was determined from pressure- 
flow characteristics using water. The mathematical 
form of the relation is a parabola, 


Ap = KF 


The performance index, P, was used to compare the 
three models and is defined as 


where F is the total flow, - is the stage efficiency, 
and 7 is the scale factor. 

Table 1 shows the performance index for each of 
the three contactors. The performance index, P, 
increases markedly as the size of the model de- 
creases. A log—log plot (not shown) indicates that P 
is proportional to 7 '*"". 
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Fig. 8 Prototype stacked-clone contactor 
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Table 1 DATA SUMMARY" 


0.5-scale unit Full-scale unit 1.4-scale unit 

Scale factor, r 0.5 1.0 1.4 
Hydraulic K 16.7 0.99 0.28 
K-r4 1.04 0.99 1.075 
Flow capacity (flooding) at 40°C and 

A/O = 3, liters/min 1.25 3.08 3.67 
Typical stage volumes 

Hydroclone cavity, underflow chamber, 

plenums, ports, and ducts, cm? 8.5 68.4 186.0 

Pump-impeller cavity and ports, cm? 41.0* 41.0 82.0* 

Pump connecting lines, cm? 50.6* 50.6 184.0* 

Total volume per stage, cm? 100.1 160.0 452.0 
Residence timej per stage at flooding, sec 1.0 3.1 { o- 
Average stage efficiency, % 85 65 74 
Performance index, P = FE/r? 4.24 2.0 1.39 


*Not scaled properly; may be varied depending upon design considerations. 
*Based on flooding data at a flow ratio (A/O) of 3 and on the volumes of the hydroclone, related 
parts, and appropriately scaled pump and lines. 


6. K. Watanabe, Studies on the Solvent Extraction of Ac- 
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NONAQUEOUS PROCESSING 


Volatility Processes 


By J. J. Barghusen 


Nitrofluor Process 


The Nitrofluor process for treating reactor fuels is 
based on dissolution of fuel and cladding in anhydrous 
HF solutions containing N,Q,, N.O3, or NOF as com- 
plexing agents. After the fuel is dissolved, the liquid 
phase is evaporated and the solid residue is contacted 
first with liquid BrF; and then with gaseous fluorine 
for recovery of uranium and plutonium, respectively. 
Recent work at Brookhaven National Laboratory 
(BNL)! has been concerned with engineering-scale 
demonstrations of the dissolution of UO, fuel clad 
with 30-mil Zircaloy in the 95°C azeotrope solvent 
(NOF -3.1HF -<0.1 N,O3) and in a NO,—HF solution. 

The dissolution tests using the 95°C azeotrope 
solvent were performed in a 16-in.-diameter Nickel 
200 dissolver using 41 gal of solvent. The experi- 
ment involved the processing of an assembly of 20 
fuel rods, 30 in. long, and containing 6.8 kg of sin- 
tered UO, pellets. The fuel bundle was supported 
within a wire basket. The dissolver solution, in this 
experiment, was the 95°C azeotrope containing an 
additional 10% free HF. The NO gas evolved during 
dissolution was vented continuously at 100 psig. 
After 20 min the solution temperature remained 
constant at 50°C and the venting of NO gas ceased, 
indicating that the Zircaloy cladding had completely 
dissolved. The solution was then slowly heated to 
60°C whereupon dissolution of the UO, pellets com- 
menced. After 2.7 hr the pellets had completely 
dissolved. 

Subsequently the liquid phase was vented and the 
solid residue, mainly NOUF,, was treated three 
times with BrF; to convert the uranium to UF¢. 
Fluorination of the salt was performed using 27-, 
13-, and 23-kg batches of BrF; at temperatures (and 
durations) of 94°C (6 hr), 127°C (13 hr), and 132°C 
(16 hr), respectively. Owing to the failure of the 
heaters on the dissolver vessel to provide a temper- 
ature of 150°C during fluorination, only one-half of 
the uranium was recovered as UF, in the treatments 
with BrF3. 
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A second experiment! involved the processing of 
16.33 kg of sintered UO, pellets, clad in Zircaloy, 
in 60 gal of 25 mole % NO,—75 mole % HF solution. 
The procedure involved the periodic addition of NO, 
to the dissolver already containing HF and the fuel 
assembly. Upon addition of half the required N,O,, 
the solution temperature and system pressure rose 
rapidly owing to reaction of the solvent with the 
cladding. Venting of the system commenced at 150 
psig; however, a peak pressure of 230 psig was 
reached before the system pressure began to de- 
crease. The remaining N,O, was added intermittently, 
and the dissolution of the fuel was finished in 7 hr. 
This experiment indicated that the dissolution of 
oxide fuel in the NO,—HF solvent canbe accomplished 
without operational problems. 

Since the development of the Nitrofluor process 
has reached the point where meaningful information 
is available for evaluating the technical and economic 
potentialities of the process, all work at BNL on 
the process has been terminated. 


Aquafluor Process 


General Electric Company has requested a con- 
struction permit and operating license for a nuclear 
fuel-reprocessing facility to be located in Grundy 
County, Ill.? This facility, Midwest Fuel Recovery 
Plant (MFRP), will utilize the Aquafluor process, 
which combines aqueous and fluoride-volatility tech- 
niques for the separation and recovery of uranium, 
plutonium, and neptunium from thermal reactor fuels 
(see Fig. 1). The plant is designed to handle 300 tons 
of irradiated uranium per year in the form of com- 
pacted UO, clad in stainless steel or zirconium alloy. 

The process employs a single-cycle solvent ex- 
traction step to separate uranium, plutonium, and 
neptunium from the bulk of the fission products. 
The actinides are stripped from the solvent. Pluto- 
nium and neptunium are recovered by anion-exchange 
techniques. The uranium effluent stream from the 
ion-exchange contactors is prepared for uranium 
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recovery by concentration of the solution in a ther- 
mosyphon reboiler and calcination of the concentrated 
uranyl nitrate hexahydrate (UNH) solution to UO. 
The calciner is an electrically heated, fluid-bed 
reactor similar to the pilot-scale unit operated at 
Weldon Spring.° Conversion of UNH to UO; is accom- 
plished at a bed temperature of 300°C using steam 
as the fluidizing gas. The UO; product from the 
calciner is mixed with inert aluminum oxide powder 
and the mixture fed by a screw conveyor to a fluid- 
bed reactor wherein the UO; is converted to UF, by 
reaction with fluorine.’ The fluorinator is 18 in. in 
diameter and is cooled by circulation of water 
through coils bonded to the wall of the reactor. 
Fluorine is fed to the reactor at a rate such that the 
off-gas stream contains about 18 vol.% fluorine, 
allowing for a normal rate of UF, production. 


The off-gas stream will contain, in addition to UF, 
and fluorine, volatile fluorides of fission-product 
elements such as niobium, ruthenium, and techne- 
tium and small quantities of volatile NpF, and PuF¢. 
Partial decontamination from these elements is ef- 
fected when the off-gas stream is passed through a 
bed of NaF pellets (2-ft bed diameter) maintained at 
400°C. Under these conditions some of the volatile- 
fluoride impurities (principally PuF,s and ruthenium 
and niobium fluorides) are sorbed on NaF, whereas 
UFs is essentially unabsorbed. The gas stream from 
the sorbent bed is cooled in a finned-tube heat ex- 
changer and then filtered by a system of sintered 
metal filters to remove particulate material. Uranium 
hexafluoride is collected in a series of manifolded 
cold traps. The traps, which are individually heated, 
are operated in a cyclic manner: on-line operation 
at low pressure, during which UF, is desublimed 
from the gas stream; and off-line operation at high 
pressure, during which the solid UF, is melted and 
drained to a product cylinder. Fluorine is recycled 
to the fluid-bed fluorinator with sufficient gas being 
bled off for disposal to maintain the desired system 
pressure. Any UF¢ in the bleed stream is recovered 
in a duplicate system of cold traps. 


Final decontamination of the UF, is effected by 
vaporizing the liquid UF, from the product cylinder 
through a 10-in.-diameter MgF, sorption bed (to 
remove NpF¢s and technetium compounds by selective 
sorption) followed by a single-stage distillation step. 
The distillation column is 6 ft long, 4 in. in diameter, 
packed with a wire-mesh material, and includes a 
reflux condenser, a heel pot, and a product con- 
denser. The system is designed to operate at 65 psi 
and 100°C. High-boiling fluoride impurities are col- 
lectec in the heel pot and are periodically removed 
and disposed of in the fluorine-waste-gas scrubber 
system (a bubble-cap tray scrubber using KOH 
solution). The purified UF, product is loaded into 
shipping cylinders using continuous weighing during 
filling of the cylinder. 
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Facilities for the production of 1500 lb of elemen- 
tal fluorine per day by fused-salt electrolysis of 
HF-—KHF, solution are provided in the MFRP. The 
fluorine-generation system employs seven 6000-amp, 
series-connected cells, each cell containing 3000 lb 
of electrolyte. The cells are water cooled to provide 
a uniform cell temperature of 83°C. 


Purification of Uranium Hexafluoride 
and Plutonium Hexafluoride 


In the fluid-bed fluoride-volatility 
reactor fuel materials, volatile UF,g and PuFs are 
produced by reaction of the fuel with gaseous fluo- 
rinating agents such as CIF, Cl1F3, BrF;, and fluorine. 
In these reactions volatile-fluoride compounds of 
certain fission products are also produced and appear 
in the UF, and PuFs, process streams. [See process 
flow sheets, Figs. 2 and 3, in Power Reactor Tech- 
nology and Reactor Fuel Processing, 10(1): 74 (Win- 
ter 1966—1967).] Some of the recent work on devel- 
oping methods for removing these compounds from 
process streams is described below. 


processing of 


REMOVAL OF IMPURITIES FROM UF; 


The use 
of solid sorbents that selectively collect fluoride 
impurities from a UFs gas stream while allowing the 
UF, to pass through the system unaffected is being 
investigated at the Oak Ridge Gaseous Diffusion 
Plant.’ A series of tests were performed to evaluate 
NaF at 344 or 400°C and MgF, at 120°C as sorbents 
for NbF;, TaF;, RuF;, SbF;, and TiF,. In these tests 
UF, gas streams containing fluorine and the fluoride 
impurity were circulated continuously through a 
packed bed of sorbent. 

Both sorbents were satisfactory in reducing the 
NbF; concentration in UF, gas streams from 460 ppm 
to less than 1 ppm prior to breakthrough. Niobium 
loadings on both sorbents averaged 5 to 7 wt.4. Some 
caking of NaF was experienced at the inlet to the 
packed bed in tests in which the NbF; concentration 
was high. This was not observed with MgF». 

Preliminary data’ suggest that MgF, at 120°C isa 
satisfactory sorbent for removing TiF, from UFs if 
the concentration of TiF, in the gas stream is less 
than 10 ppm. Titanium concentrations in the gas 
leaving the MgF, trap were less than 1 ppm. Mean- 
ingful data have not yet been obtained for UF, streams 
containing higher concentrations of TiF4. 

Sodium fluoride at 344°C appears to be a good 
sorbent for RuF;, reducing the ruthenium concen- 
tration in the gas stream from 250 to 2 ppm. The 
volatile ruthenium fluoride used in this study (pre- 
sumably RuF;) was prepared by direct fluorination 
of the metal, and therefore the effectiveness of NaF 
for trapping ruthenium can not be extended with 
certainty to other volatile ruthenium fluorides or 


Evaluation of NaF and MgF> as Sorbents. 
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oxyfluorides produced by fluorination of oxides. It 
was observed that the NaF—RuF; complex was un- 
stable and tended to decompose and desorb upon 
exposure to air. This behavior precluded accurate 
determination of ruthenium loading and identification 
of the complex.* 

Effective trapping of TaF; (<1 ppm TaF; in the exit 
gas) was achieved‘ by the use of MgF, at 120°C or 
NaF at 344°C. Preliminary results indicate that 
MgF, at 120°C is a satisfactory sorbent for removing 
SbF; from UFs. 

In general, the quantity of uranium picked up by 
the sorbent in any of the tests was less than 2 x 107° 
wt.% of the total quantity of uranium passed through 
the sorbent bed. 


Revised Specifications for UF;. Uranium hexaflu- 
oride derived from the reprocessing of highly irra- 
diated power-reactor fuels will contain certain im- 
purities at levels not present in UF, derived from 
both normal uranium and low-burnup reactor fuels. 
These impurities are 23217 and alpha-active trans- 
uranium elements. To establish whether the 1958 
specification limits for these materials were neces- 
sary, AEC appointed a Uranium Hexafluoride Specifi- 
cations Study Committee in January 1966 to reexamine 
all UF, specifications with respect to the need for 
these specifications.®° The examination indicated that 
all specifications were required for safe and econom- 
ical operation of the diffusion plants and that all 
specification limits were economically achievable by 
manufacturers of normal UFs,. However, the study 
showed that the limit for 7°*U of 11 parts per billion 
parts *“*u and the transuranic alpha-activity limit 
of 150 dis/min per gram of uranium, which were 
previously in effect, could not be met in UF, from 
spent power-reactor fuels. On the basis of the rec- 
ommendations from the committee, AEC announced 
a 10-fold relaxation in the *U and transuranic 
alpha-activity limits.° The specification limits for 
UF, delivered to AEC, effective July 1, 1967, are 
presented in Table 1. 

For several years, UF, containing very low con- 
centrations of neptunium and plutonium was fed to 
the Paducah Gaseous Diffusion Plant. Both neptunium 
and plutonium have been found on plant-equipment 
surfaces but neither has been detected in product 
streams. Since equipment contamination is a cumu- 
lative problem, the level of contamination will in- 
crease as the total quantity of neptunium and pluto- 
nium fed in the UF, increases; thus dilution by 
uncontaminated UF, is not a factor. Assessment of 
methods for removing these species from UF, awaits 
the development of accurate sampling techniques for 
trace amounts of the materials.° Preliminary tests 
have shown that alkaline-earth and alkali-metal flu- 
orides are promising agents for removal of NpF, 
and PuF, from UF,. Table 2 lists the results of PuF, 
sorption on metal fluorides at several temperature 


Table 1 SPECIFICATIONS FOR URANIUM HEXAFLUORIDE 
DELIVERED TO THE AEC® 


(Effective Date July 1, 1967) 








Numerical 
Item* value 

Maximum vapor pressure of filled container at 200°F 

in pounds per square inch, absolute 75 
Minimum weight percent of uranium hexafluoride in 

material 99.5 
Maximum mole percent of hydrocarbons, chlorocar- 

bons, and partially substituted halohydrocarbons 0.01 
Maximum number of parts of elements indicated per 

million parts of total uranium 

Sb 1 

Br 5 

Cl 100 

Nb 1 

Pp 50 

Ru 1 

Si 100 

Ta 1 

Ti 1 

otal of elements forming nonvolatile fluorides 

(having a vapor pressure of 1 atm or less at 300°C), 

e.g., Al, Ba, Bi, Cd, Ca, Cr, Cu, Fe, Pb, Li, Mg, Mn, 

Ni, K, Ag, Na, Sr, Th, Sn, Zn, and Zr 300 
Maximum number of parts of elements or isotopes 

indicated per million parts of 25U 

Cr 1500 

Mo 200 

WwW 200 

V 200 

233] 500 

2s2y 0.110 
Maximum thermal neutron absorption of total im- 

purity elements as equivalent parts of boron per 

million parts of total uranium 8 
Maximum total of gamma activity due to fission 

products and 2°7U as percent of gamma activity 

of aged natural uranium and as measured in a 

high-pressure ionization chamber (Drawing D- 

AWM-8796 of Nuclear Division, Union Carbide 

Corporation) 20 
Maximum beta activity due to fission products 

as percent of beta activity of aged natural uranium 10 
Maximum number of parts of plutonium per billion 

parts of total uranium 1 


Maximum alpha activity from all transuranic elements 
in disintegrations per minute per gram of total 
uranium 1500 





*All specification analyses on uranium hexafluoride shall be 
performed on samples removed in the liquid state from each cyl- 
inder while its contents are liquid and homogeneous 


levels.° Although these tests involved the sorption of 
PuF, from a flowing fluorine stream at higher plu- 
tonium concentrations than will be encountered ina 
diffusion-plant feed, they indicate that several sor- 
bents merit further studies. 


NbFs and RuFs Removal by Metal Fluorides. At 
the 36th International Congress on Industrial Chem- 
istry, Prusakov’ reported the results of a study on 
the removal of NbF; and RuF; from UF, gas streams 
by sorption on metal fluorides. The most active 
sorbents for NbF; were NaF at 300°C, LiF at 100°C, 
and BaF, at 100 to 200°C. The interaction of NbF; 
with NaF at 100, 200, and 300°C produced sodium 
heptafluoroniobate (Na,NbF7), whereas hexafluoroni- 
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Table 2. TEMPERATURE DEPENDENCE FOR RETENTION 
OF PuFg ON GROUP IA AND ITA METAL F LUORIDES® 





Percentage of plutonium retained at temperature, °C 


Sorbent 100 200 300 400 500 600 700 








LiF 29.4 46.8 87.5 88.7 96.9 89.7 90.4 
Nak 98.6 99.8 99.6 99.92 99.94 99.97 98.2 
KF 97.1 98.5 97.6 97.5 94.2 94.4 98.0 
RbF 83.6 90.4 93.8 99.9 99.9 98.5 
CsF 99.6 98.6 97.9 99.8 95.8 
Bek, 29.1 20.1 76.6 70.9 71.0 73.4 87.1 
MgF, 98.8 98.2 96.3 96.9 98.6 98.6 94.5 
CaF, 99.99 99.995 99.994 99.997 99.98 99.996 
SrFy 99.95 99.97 99.7 98.9 99.97 99.96 99.97 
BaF, 99.4 99.95 99.95 99,99 99.97 93.93 99.99 


obates were produced by reaction of NbF; with LiF, 
BaF., MgF., and CaF». The reaction of RuF; with 
NaF, LiF, and BaF, at 100 to 300°C and with MgF, 
and CaF, at 100°C produced hexafluororuthenates of 
the corresponding metals. At 300°C, RuF; may be 
desorbed from the sodium fluoride complex in a 
fluorine gas stream. The dissociation pressures of 
several complex fluorides of niobium and ruthenium 
were measured by transpiration techniques, and the 
data are presented in Table 3 as constants for the 
equation log P = A—B/T, where P is the equilibrium 
pressure in torrs of the volatile fluoride over the 
solid phase, and J is °K. 


Table 3) THERMAL -DISSOCIATION PRESSURE OF 
fF LUORONIOBATES AND FLUORORUTHENATES 
AS A FUNCTION OF TEMPERATURE? 
[log P (torr) = A—B/T (°K)] 
Temperature 

System interval, °C A B 
LiF —NbFs 210-290 5228 2456 
NaF —NbF; 330-480 4410 2437 
MgF,—NbF; 120-218 4881 1659 
CaF,—NbFs 87 — 223 3039 864 
BaF, —NbF; 88 —247 1995 464 
LiF —RuF; 229-360 5992 2878 
NaF —RuFs 313 —490 7015 4004 
MgF, —Rufk’s 165-301 3821 1445 
CaF, —RuFs 196 —375 2265 845 
BaF, —RuFs 238-415 3994 2312 


METHODS FOR PURIFICATION OF PuF; 


In the fluoride-volatility process flow sheets based 
on the use of interhalogens for selective fluorination 
of UF,, plutonium is recovered as PuF, by fluorina- 
tion with fluorine. Several fission products that form 
volatile fluorides will accompany the PuF¢,, namely, 
ruthenium, antimony, and niobium; and the PuF¢, 
stream must be processed further for the prepara- 
tion of a purified plutonium product. 


Sorption of PuFs on LiF. Workers at Oak Ridge 
National Laboratory’ are investigating the use of LiF 
as a selective sorbent for PuFs since this fluoride, 
unlike NaF and other metal fluorides, permits nearly 
quantitative desorption of PuFs at elevated tempera- 
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tures. The rate of PuF, sorption on LiF was observed 
to be affected by the surface area of the sorbent. 
High PuF, sorption was achieved at 175°C with LiF 
having a surface area >20 m’/g, whereas a tempera- 
ture of 300°C was required to achieve a similar 
sorption rate with low-surface-area (3 m?/g) LiF. 
The complex formed at 100°C with high-surface-area 
LiF has an average plutonium valence of 4.55. Tests 
to demonstrate the effectiveness of LiF in recovering 
PuF, from simulated process gas streams are being 
carried out in fluid-bed systems. 


Sorption of PuFs on NaF, Experimental results re- 
ported by Prusakov’ and by Argonne National Lab- 
oratory (ANL)® have indicated that the reaction 
between PuFs, and NaF at 100 to 150°C produces the 
complex compound 3NaF- PuFy. Other experiments at 
ANL? have shown that the complex salt 7NaF -6PuF,, 
prepared by aqueous chemistry, upon contact with 
fluorine at temperatures up to 550°C releases PuF, 
to produce the compound 2NaF-PuFy,. This material 
is highly stable in fluorine up to 550°C, and further 
recovery of plutonium from this compound by vola- 
tility techniques would be difficult. 

A technique for purifying PuFs, process gas streams 
based on the sorption of PuF, on NaF at 100°C is 
being explored at ANL.!’ Under the sorption condi- 
tions considered for PuFs, UF,, NpF., Tc, and small 
amounts of Ru, Nb, and Sb are also sorbed. Removal 
of uranium, neptunium, and technetium from the NaF 
bed is achieved by passing fluorine through the bed 
at 400°C. The resulting NaF bed, containing pluto- 
nium, is leached with anhydrous HF at 20°C. Sodium 
fluoride is soluble in liquid HF. Plutonium tetrafluo- 
ride, which is insoluble, appears as a precipitate 
and is recovered by filtration of the solution. 


BrF; Reduction of PuFs. An alternative method!” 
for purifying PuF, streams involves the selective 
reduction of PuF, to PuF, with BrF3. The procedure 
involves contacting the mixture of PuF¢., actinide 
hexafluorides, and fission-product fluorides with BrF3 
for 1 hr at 80°C. The reaction vessel is then evac- 
uated and the temperature of the vessel increased 
to 150°C to volatilize the fission products and bromine 
fluorides. Remaining traces of uranium, neptunium, 
and fission products are removed from the PuF, 
product by contact with BrF; gas at 300°C. 

Experiments have been performed!','? to evaluate 
the reduction of PuF, from mixtures containing gram 
quantities of UF,, NpF,., SbF;, NbFs, and RuF;. Re- 
duction of PuF, to PuF, was complete. The major 
portion of uranium, antimony, and niobium was sep- 
arated from plutonium in the evacuation step. Even 
higher decontamination from these elements was 
achieved in the step in which the PuF, was contacted 
with BrF; for 6 hr at 300 to 400°C. 

Essentially no decontamination from ruthenium 
was observed after the reduction-evacuation step. 
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The decontamination factor for ruthenium after the 
BrF; fluorination step was approximately 540. No 
significant separation of neptunium was achieved. 
Preliminary experiments!” using tracer quantities of 
niobium and ruthenium showed unsatisfactory removal 
of these elements from plutonium. 


Thermal Decomposition of PuFs. Thermal decom- 
position of PuF, to PuF, at 300°C is being considered 
as a technique for separating plutonium from mix- 
tures containing UF, and other fluorides. Since nep- 
tunium may be present in these mixtures as NpF¢, 
the behavior of NpF, during thermal decomposition 
of PuF, was determined at ANL.'* Gaseous mixtures 
of NpF,s and PuFs, were circulated through a vessel 
at 300°C. It was observed that a small fraction of 
the neptunium (~0.2 wt.%) codeposited with the plu- 
tonium. The concentration of neptunium in the de- 
posited solids was constant throughout the entire 
experiment. The fraction of the total neptunium in 
the gas mixture that codeposited with plutonium was 
found to be dependent upon the total quantity of plu- 
tonium deposited. 

Thermal decomposition of PuF, has been demon- 
strated at ANL on an engineering scale using a fluid- 
bed reactor to carry out the decomposition reaction. 
Three experiments’ were made to examine the 
feasibility of separating plutonium from UF,—PuF, 
mixtures by thermal decomposition. In the process 
a mixture of gases is passed through a fluidized bed 
of sintered, refractory alumina (a bed of PuF, might 
be used in a large-scale continuous operation) at 
300°C, and PuF, deposits on the bed material. In 
each experiment 10 kg of a PuF,—UF, mixture (up 
to 0.4 wt.% plutonium), diluted to 40 vol hexafluo- 
ride with nitrogen, was processed in 8 hr. 

The results of the experiments indicated that the 
decomposition process is very efficient as evidenced 
by the low plutonium content of the UF, off-gas 
stream. On the basis of analyses of process streams, 
a plutonium separation of about 99.9% was achieved. 
The final alumina bed contained 1.46 wt.% plutonium 
and only 0.19 wt.% uranium. The decomposition of 
PuF, to PuF, occurred preferentially on the sintered 
alumina-bed particles rather than in the gas phase. 

Since the UF,—PuF, feed stream to the fluid-bed 
decomposer was derived from fluorination of sintered 
UO,—PuOy, pellets containing synthetic fission prod- 
ucts, the mixed feed contained some ruthenium 
fluoride. Analyses of the bed material and process 
gas streams indicated that no ruthenium was depos- 
ited with the PuF, on the bed and that the ruthenium 


passed through the decomposer with the UF, (Ref. 
15). The ruthenium content of the PuF,—alumina bed 
was 0.2 ppm. 
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NONAQUEOUS PROCESSING 


Compact Pyrochemical Processes 


By W. E. Miller and R. K. Steunenberg 


Pyrochemical processes employ high-temperature 
nonaqueous methods for separations of fissile, fertile, 
and fission-product elements in irradiated reactor 
fuels. In most cases liquid metals and/or salts are 
used as process solvents. The required separations 
may be effected by selective volatilization, precipita- 
tion from a solvent, liquid metal—salt extraction, or 
electrolysis. Processes for ceramic (oxide or car- 
bide), metallic, and molten-salt fuels are currently 
being developed. Much of the development work is 
directed toward the reprocessing of fast breeder re- 
actor fuels. Processes of this type are expected to 
be capable of handling high-burnup, short-cooled 
fuels in compact equipment. 


Salt-Transport Process 
for Fast Breeder Reactor Fuels 


Current work at Argonne National Laboratory is 
directed toward the development of a pyrochemical 
process for the recovery of irradiated core and 
blanket fuel from fast breeder reactors. A conceptual 
flow sheet for this process was presented in detail in 
an earlier issue of Power Reaclor Technology and 
Reactor Fuel Processing.' The process represented 
by the flow sheet is for UO,—PuO, fuels clad with 
stainless steel. The oxide is first declad by selec- 
tively dissolving the stainless steel in liquid zinc at 
800°C. The oxide is then separated from the zinc 
phase, suspended in a molten salt (47.5 mole % 
MgCl,—47.5 mole % CaCl,—5.0 mole % CaF,), and 
reduced by contacting the salt with a liquid Cu— 
33 wt.% Mg alloy at 800°C. The gaseous fission 
products are released to the system atmosphere dur- 
ing the reduction step. The alkali, alkaline-earth, and 
rare-earth fission products and the MgO by-product 
of tne reduction reaction appear in the salt, which is 
discarded as waste. The metallic uranium, plutonium, 
and the remaining noble and refractory fission prod- 
ucts appear in the liquid-metal phase. The plutonium 
is then separated from the uranium and the fission 


products by cycling a new molten salt (50 mole % 
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MgCl,—30 mole % NaCl—20 mole % KCl) between the 
Cu—33 wt.% Mg alloy and a Zn—2 wt.% Mg alloy at 
600°C. With this salt-transport step, the plutonium 
can be selectively transferred to the Zn—2 wt.% Mg 
alloy. The zinc and magnesium are removed from the 
metallic plutonium product by vacuum distillation (re- 
torting). The uranium, which is present in the Cu— 
33 wt.% Mg alloy as precipitated metal. is recovered 
by decanting the supernatant liquid-metal alloy and 
washing the precipitate with liquid magnesium. Re- 
sidual magnesium is removed from the uranium by 
retorting. Mixed UO,—PuO, is prepared from the 
product metals by direct oxidation. 


DEMONSTRATION OF SALT-TRANSPORT PROCESS STEPS 


Recently reported work on the salt-transport pro- 
cess?.3 has been concerned primarily with a demon- 
stration of the flow sheet, with the exception of the 
decladding and oxide reconversion steps. The equip- 
ment‘ used in this demonstration, which consists of 
two resistance-heated reaction vessels interconnected 
by a heat-transfer line to permit the movement of 
process solvents between the vessels, is shown in 
Fig. 1. The starting material consisted of 233 g of 
PuO,, 458 g of UO,, and 52 g of oxides of inactive 
fission-product elements. The oxide-reduction and 
salt-transport steps of the process were easily done. 
In the salt-transport step, which was described in a 
previous issue of Reactor and Fuel-Processing Tech- 
nology,® preliminary analyses indicate that 99.4% of 
the plutonium was transported in 14 cycles. An addi- 
tional processing step was carried out to remove any 
copper that may have been transported with the plu- 
tonium product. In this step ZnCl, was added to the 
Zn—2 wt.% Mg alloy to convert the plutonium to 
PuCl;, which was extracted into the salt phase. The 
plutonium was then extracted from the salt by a Zn- 
4 wt.% Mg alloy. The plutonium recovery in this 
step was 99.9%, and the copper content of the Zn— 
4 wt.% Mg alloy was 20 ppm, which corresponds to 
the amount of copper in the reagents making up the 
alloy. 
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Equipment for testing salt-transport process.* 


Fig. 1 


TWO-PHASE THERMAL-CONVECTION LOOP 
FOR CORROSION TESTING 


A two-phase thermal-convection loop made from 
1-in.-diameter tubing is being developed for corro- 
sion tests of potential container materials for pyro- 
chemical processes.” The equipment consists of two 
loops, one on top of the other, in the form of a modi- 
fied figure “8.” The upper loop is filled with salt, and 
the lower one is filled with metal. The interface 
between these two immiscible phases is contained in 
a horizontal leg common to both loops. Temperature 
gradients are maintained between the vertical legs of 
each loop so that countercurrent thermal-convection 
flow occurs between the metal and salt phases. 

In the first test the loop was constructed of type 304 
stainless steel. The purpose of this test was to study 
the operating characteristics of the double-loop de- 
sign. The metal phase was Cu—33 wt.% Mg, and the 
salt phase was 50 mole % MgCl,—30 mole % NaCl-— 
20 mole % KCl. The salt-phase temperatures were 
819°C in the hot leg and 736°C in the cold leg. The 
metal-phase temperatures were 830°C in the hot leg 
and 650°C in the cold leg. 

After 7 days of operation, the circulation rate of 
the metal phase decreased. The loop was then shut 
down and sectioned, and deposits of iron and chro- 
mium from the wall of the loop were found in the cold 
leg of the metal loop. This mass transfer of wall 
material from the hot leg to the cold leg was similar 
to that observed when liquid metals, such as bismuth, 
lead, and mercury, are contained in chrome-—nickel 
steels with a temperature gradient. 

This test has demonstrated the feasibility and 
usefulness of the double-loop design for studying 
corrosion in two-phase systems. Loops made from 
more promising container materials, suchas tantalum 
and niobium, will be used to study corrosion of these 
materials. 





NEPTUNIUM SOLUBILITY STUDIES 


Experiments were performed to determine the solu- 
bility of neptunium in liquid cadmium and the distri- 
bution coefficient of neptunium between a molten salt 
(50 mole % MgCl,—30 mole % NaCl—20 mole % KCl) 
and various cadmium— magnesium alloys.’:3 With this 
information the possibility of recovering neptunium in 
pyrochemical processes for fast breeder reactor 
fuels was investigated. Initial results of the solubility 
experiments indicated a solubility of about 0.3 at.% 
neptunium in liquid cadmium at 393°C. The values of 
the distribution coefficient* at 650°C ranged from 3 
to 0.15 for Cd—Mg alloy compositions ranging from 
0.7 to 81 wt.% magnesium, respectively. These re- 
sults indicate that a neptunium recovery step may be 
feasible. 


Processes for Molten-Salt Breeder 
Reactor (MSBR) Fuel 


Oak Ridge National Laboratory (ORNL) is develop- 
ing a thermal breeder reactor in which the core fuel 
is a molten fluoride salt (66 mole 7% LiF-—34 mole % 
BeF,) containing *°U or *U fluoride as the fissile 
material and ZrF, as the stabilizing agent to prevent 
uranium oxide precipitation. The blanket, which con- 
tains fertile thorium, is a molten fluoride salt of the 
nominal composition 71 mole % LiF—2 mole % BeF,-— 
27 mole % ThF,. Nonaqueous techniques are being 
developed to reprocess both the core and the blanket. 
The objective of reprocessing the core fuel is to re- 
move rare-earth fission products, which are thermal- 
neutron poisons. The objective of blanket processing 
is to recover fissile 7*°U and bred ***Pa, which decays 
to °3U, Research and development work on these pro- 
cesses is discussed in a recent ORNL report.°® 

The major steps® in the present conceptual process 
for treating molten-salt breeder reactor fuel are 
shown in Fig. 2. Core fuel is first contacted with 
fluorine to oxidize uranium to UFs, which volatilizes 
from the liquid salt. Fluorination also removes certain 
fission-product fluorides, such as fluorides of Mo, 
Te, Ru, Zr, and Nb, which volatilize along with the 
UF;. The UF, stream is passed through salt-sorber 
beds to remove fission-product fluorides from the 
UF,. The barren fuel-carrier salt leaving the fluori- 
nator is LiF—BeF, containing rare-earth fission 
products. The carrier salt is purified by vacuum 
distillation of LiF—BeF»,. The still residue consists 
of a highly concentrated solution of rare earths in 
LiF, which is discarded as waste. The core fuel is 
reconstituted in a single operation in which UFs and 
H, are introduced into the purified molten salt and 
the UF, is reduced to UFy. After some feed adjust- 





«Distribution coefficient (K,) = (wt.% neptunium in salt 
phase)/(wt.% neptunium in metal phase). 
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ment, the processed fuel salt is returned to the re- 
actor core. 

The feasibility of distilling fluoride salts was 
established in batch laboratory experiments by Kelly.’ 
Equipment has now been designed for an engineering- 
scale demonstration of vacuum distillation of molten- 
salt reactor fuel. This demonstration will show that 
molten salt containing fission products can be fed 
continuously to the still at the same rate at which it 
is being distilled with the simultaneous accumulation 
of fission products in the bottoms. The still can also 
be operated batchwise to concentrate fission products 
in some small fraction of the original charge. The 
vacuum-distillation equipment is shown in Fig. 3. 
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Fig. 3 Vacuum distillation® of LiF ~-BeF>—ZrF4. 


To design a still requires knowledge of the relative 
volatilities of the rare-earth fluorides and LiF, which 
is the major constituent in the still pot. The volatili- 
ties of four of the rare-earth trifluorides, which have 
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been obtained from recent experiments, relative to 
that of LiF at 1000°C and 0.50 mm Hg are listed as 
follows: 


Rare-earth Mole Relative volatility 


fluoride fraction with respect to LiF 
CeF 0.02 Sx a0e 

LaF, 0.02 320 

NdF 0.05 6 x 107+ 
SmF 0.05 2x10 





These relative volatilities show that rare-earth fluo- 
ride removals of the required efficiencies can be 
obtained in a still of simple design without rectifica- 
tion. 

Two other methods of separating rare-earth fission 
products from the molten fuel-stream carrier salt 
are under investigation: (1) the rare-earth fission 
products may be reduced and extracted into a molten- 
metal phase that is immiscible with the salt or (2) they 
may be removed by the formation of insoluble solid 
solutions with UF3. 

When the molten fuel-stream carrier salt is con- 
tacted with a Bi—Li alloy, rare-earth fluorides are 
reduced, and the resulting rare-earth metals deposit 
in the alloy. Experimental results have indicated that 
much more than the theoretical amount of lithium 
disappears from the metal alloy when the alloy and 
the salt are contacted. It is suspected that lithium 
reduces some of the BeF, in the salt. 

With the other method of rare-earth removal, the 
rare-earth fluorides can be precipitated with UF; as 
solid solutions by adding solid UF; to the molten salt. 
Results obtained for the removal of NdF3 are illus- 
trated in Fig. 4. Material-balance calculations for 
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Fig. 4 Removal’ of NdF3 from LiF—BeF» (66 to 34 mole %) 
at 550°C by precipitation on solid UF}. 
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the system indicate that the composition of the solid 
phase, resulting from the addition of UF3, varied 
from 0.54 to 0.22 mole % NdF; in UF 3. Subsequent 
changes in the concentrations of NdF; and UF; in the 
salt mixture, as the temperature was varied from 
550 to 800°C, indicated a constant solid-phase com- 
position of 0.22 mole % NdF; in UF3. 


Measurements of the solubilities of rare-earth 
fluorides in mixtures of molten LiF—BeF, have been 
resumed to supplement earlier data. The purpose of 
this experimental program is to investigate the 
behavior of those rare earths which are of interest to 
fuel-reprocessing studies for the reference design 
MSBR and whose solubilities in the proposed MSBR 
fuel solvent have not been previously measured. The 
solubility of SmF3; in 66 mole % LiF—34 mole % BeF, 
was determined for comparison with previous SmF3 
solubility data in a similar fluoride solvent (62.8 
mole % LiF-—36.4 mole % BeF,—0.8 mole % UF4,) and 
with anticipated measurements of SmF, solubilities 
in the LiF—BeF, solvent. Experimental values for the 
solubility of SmF, at 600, 700, and 800°C were 0.013, 
0.024, and 0.040 mole fraction, respectively. Solubili- 
ties of NdF; in the same solvent were 0.010, 0.019, 
and 0.035 mole fraction at 600, 700, and 800°C, re- 
spectively. Heats of solution calculated from the 
linear dependence of the logarithm of the solubility 
values on the reciprocal of the absolute temperature 
were 10.5 kcal/mole for SmF3 and 11.5 kcal/mole for 
NdFs3. 

The solubility values for SmF3 in LiF—BeF, coin- 
cide only at 650°C with those obtained previously for 
the LiF—BeF,—UF, solvent. The heat of solution of 
SmF3, based upon the previous data, is about 13.9 
kcal/mole. Although no discrepancies in the two ex- 
perimental procedures are apparent, further studies 
would be required to investigate possible solvent 
effects. 

Thorium is the fertile material in the MSBR blanket. 
Fissile *°U is produced in the blanket by the reac- 
tions 


929 5 rd ¢ a 
2327 +  — 23H 233g B 23317 


The concentration of 7**U in the blanket must be main- 
tained at a low level to avoid appreciable parasitic 
neutron captures. A recent issue of Reactorand Fuel- 
Processing Technology® reviewed processing methods 
for removal of protactinium from blanket salt. 

In one of the processes under study, blanket salt is 
contacted with Th— Bi amalgam to reduce the protac- 
tinium fluoride to metal. In another process solid 
thorium metal is added to the blanket salt, and the 
protactinium metal formed by the reduction of the 
fluoride is collected on steel wool. Difficulties in 
collecting the reduced protactinium metal have been 
experienced in both processes. Work on these pro- 
cesses is continuing. 


Miscellaneous Developments 


Zone-melting experiments’ were conducted at Kern- 
forschungszentrum, Karlsruhe, Germany, to inves- 
tigate the removal of fission products and other 
impurities from UCl,. In all cases the impurities 
distributed strongly toward the molten zone, thus 
moving toward the end of the sample. The separation 
factors (C/C,) were between 10-4 and 5 x 10-2 and 
were in the following order: RuCl; < BaCl, < CsCl < 
CeCl; < NaCl < ZrCly. (The separation factor in this 
case is defined as the ratio of concentration of solute 
after zone melting to the concentration of solute prior 
to zone melting.) The values were computed for condi- 
tions after 20 passes of the molten zone at 10 mm/hr. 
Separations with tracer concentrations showed the 
same results as those obtained with higher concentra- 
tions (0.1—1 wt.%). Zone-refining experiments with a 
UCl,—KCl eutectic yielded less favorable results. 
Separation factors under comparable conditions were 
less favorable by one to three orders of magnitude. 
For instance, for the solute RuCl;, C/C, values were 
2.0 x 10-4 and 0.12 in UCly and UCl,—KCl, respec- 
tively. In the UCl,—KCl eutectic, zirconium tetra- 
chloride was the only compound that distributed 
toward the solid phase, and it did so only when it was 
present in macro concentrations. 

In other studies at Karlsruhe,’ zone melting and 
column crystallization were compared as separative 
methods in the system KNO;—NaNO3, the binary 
eutectic KCl— LiCl, the ternary eutectic MgCl,—NaCl— 
KCl, and the binary eutectic UCl,—KCl. The results 
showed that both purification methods are, in prin- 
ciple, satisfactory for inorganic fused-salt separa- 
tions. Zone melting appears to be well suited for the 
purifications of certain substances, but it requires a 
long period of time and only small amounts of mate- 
rial can be handled. Column crystallization, on the 
other hand, is rapid and is suitable for large through- 
puts. It offers the best separation efficiency in those 
systems in which the components are present in 
large concentrations. However, serious mechanical 
problems are encountered in applying this method to 
melts that have high freezing points or are highly 
corrosive. 

Glasses, which have been used for cation exchange 
in molten-salt solvents, show a remarkable selec- 
tivity in certain systems. Investigations were con- 
ducted on the distribution behavior of radioactive 
tracers between molten salts and clean crushed 
glasses.'! In a system in which the glass composition 
was 21.8 mole % Na,O—60.2 mole % B,O;—18.0 
mole % NaCl and the salt was molten NaNOs, order- 
of-magnitude values of K,* at 550°C for the solutes 
Cst, Ba’*, and Eu3+ were 107°, 1, and 10°, respec- 
tively. 


0 





*Kq is defined as (specific activity in glass phase)/(spe- 
cific activity in salt phase) for the adsorption process. 
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WASTE DISPOSAL 


Progress in Waste-Disposal Research 


and Development 


By Phillip Fineman 


Conversion of High-Level-Activity 
Wastes to Solids 


FLUID-BED CALCINATION 


Radioactive aqueous wastes resulting from the repro- 
cessing of spent reactor fuel are being converted to 
a granular, free-flowing solid in a heated, fluid-bed 
calciner in the Waste Calcining Facility (WCF) of the 
Idaho Chemical Processing Plant (ICPP).' The cal- 
cined solids are removed continuously from the cal- 
ciners and transported pneumatically to underground 
storage bins. Reductions in volume (iiquid to solid) 
of 9 to 11 are being achieved.!:* Presented in the 
following paragraphs are further results ofthe second 
radioactive-waste-processing campaign, currently in 
progress, and a discussion ofthe operating difficulties 
encountered3—!3 Also discussed are possible modifi- 
cations to the WCF which would improve its operating 
reliability and increase throughput rates. In addition, 
laboratory and pilot-plant investigations related to 
the following are described: calcination of zirconium 
fluoride wastes, storage of calcine, high-temperature 
performance of various stainless steels in NaK, and 
a valveless control system. 


Second Waste-Processing Campaign. Since the 
last report on the current processing campaign, ' the 
WCF has continued to calcine successfully the aqueous 
plant wastes presently being produced and stored at 
the ICPP.°- The wastes calcined during this period 
were aluminum nitrate solutions containing high con- 
centrations of nitric acid and acid-deficient aluminum 
nitrate solutions containing high concentrations of 
ammonium nitrate. Some 372,000 gal of these two 
wastes were converted into 4800 cu ft of solids. In 
general, operation of the WCF was highly satisfac- 
tory. Operating problems, which are discussed below, 
appear to be readily solvable. Preconcentrating the 
blended wastes (to 85% of the original volume) made 
it possible to increase the net throughput of the 
calciner by about 15% (based on the rate of waste 


depletion from the tank farm)?’"!-!2, The overall net 


feed rate* has averaged about 70 gal/hr, which ex- 
ceeds the design feed rate of 60 gal/hr. The off-gas 
from the calciner was effectively decontaminated; the 
release rates to the atmosphere of "Sr and !8Ru re- 
mained consistently well below the established safe 
release limits. 


Operating Problems. Although a greater number 
of operating problems were encountered in this cam- 
paign’ than in the first one,” the relative ease with 
which these problems can be solved has demonstrated 
the considerable flexibility of the WCF. The principal 
problems encountered and the means for their solu- 
tions are as follows: 

1. Plugging of the Feed System. Sand, present in 
the alumina starting-bed material, and silica gel col- 
lected in the feed evaporator and then released to the 
feed stream, caused recurring feed-control problems 
early in the present campaign.° During the first few 
days after initial start-up, finely divided sand was 
elutriated from the fluidized bed and carried into the 
WCF scrubbing system. This sand then found its way 
into the feed system by way of the scrub solution, 
which is recycled to the calciner. Here the sand 
tended to stop or impede the flow of the feed solu- 
tions. In the future it is planned to use bed sub- 
stances other than sand to eliminate this problem. 
One of the substances being considered is dolomite, f 
which is soluble in the scrubbing solution. The 
results of recent pilot-plant work***» have shown 
that dolomite can be used as a starting-bed material. 
Pilot-plant work is also under way to identify the best 





*The feed rate for the WCF calciner is also stated in 
terms of gross values. The gross feed rate has averaged 
about 93 gal/hr; it includes about 5% jet dilution that re- 
sults during transfer of the waste solutions to the WCF 
and about 22.5 gal/hr of recycle from the off-gas scrubbing 
system.’ 

*Dolomite is composed of calcium and magnesium car- 
bonates. 
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means of modifying the WCF feed system to eliminate 
plugging tendencies when feed solutions contain signif- 
icant quantities of solids. 

2. Difficulties in Processing of Al(NO3)3—NH,NO3 
Wasles 

Early in July 1966, while processing aluminium nitrate 
waste containing a high concentration of ammonium ion, 
the pressure drop across the final filters increased about 
4 inches, to 10.5 inches of water, during a four-day 
period of operation. The off-gas filters were replaced, 
and ammonium-free waste was started into the WCF be- 
cause it was believed that the rapid rise in pressure drop 
was related to the ammonium waste. However, a rapid 
pressure drop rise occurred with the new filters until the 
temperature of the off-gas entering the filters was de- 
creased from 85 to 80°C. After the temperature change, 
a dramatic four-inch decrease in pressure drop across 
the filters was observed. An analysis of the operation 
indicates that trace quantities of ammonium nitrate (first 
from the waste during calcination and later —-while am- 
monium-free waste was being calcined—from residual 
material in the calciner, scrubbing system, and/or silica 
gel absorbers in the off-gas train) may deposit on the 
filters and form crystals to partially plug the filter gas 
passages. The connection between temperature and pres- 
sure drop could be due to a temperature effect on the 
size of ammonium nitrate crystals. At 84.5°C, a marked 
change occurs in dimensions of ammonium nitrate 
erystals.¢) 

3. NaK Heating System Problems. During the 
first and second processing campaigns, problems in- 
volving the NaK heating system have caused the 
greatest amount of downtime.° Several weld failures 
on a section of piping near one of two expansion tanks! 
and a failure in one of the duplextubes in the calciner 

. : 5 
vessel have necessitated shutdowns for repairs.” In 
all cases a thorough review of the circumstances has 
linked the failures to stress conditions related to 
the design of the system. Some of the sources of ex- 
cessive stress were eliminated by removing one 
expansion tank and making some piping changes. A 
long-range study is now under way to determine the 
best design to replace the existing heating system. 
Operating experience has shown thatthe entire system 
can be simplified and that greater reliability and in- 
creased capacity can be designed into a new system 
without any sacrifice in safety. 


Possible Modifications tothe WCF. The following 
modifications for improving the WCF operating reli- 
ability and processing flexibility are under consider- 
ation:!! 

1. NaK Heating System. Several methods of heating 
the WCF are being examined for the purpose of finding 
a replacement for the liquid-metal heating system. 
These methods include electric heaters, directly 
fired tubes using combustion gases, fluidizing air at 
very high temperatures, and in-bed combustion of 
hydrocarbons. The most promising method, which is 
currently under investigation, is in-bed combustion 
of kerosene. A 10-day run’ made in the 12-in.- 
diameter pilot-plant calciner showed that in-bed 
burning provides a satisfactory means of heating 
while calcining either aluminum nitrate or zirconium 
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fluoride wastes. Equipment inspection after the run’ 
revealed no obvious corrosion effects. Analysis of 
run results is under way. 

2. Reduction of Solids Carry-over. Essentially all 
solids elutriated from the calciner bed are removed 
from the off-gas in the calciner cyclone or in the 
scrubbing system.!! About 20% of the elutriated solids 
enter the scrubbing system where they dissolve and 
periodically are recycled to the feed, thus reducing 
the overall calciner throughput. Elimination of the 
solids from the off-gas by dry filtration would greatly 
reduce, if not eliminate, the recycle of scrubbing so- 
lution. Currently, stainless-steel filters are being 
considered as replacement for the calciner cyclone. 
Corrosion tests on applicable sintered stainless- 
steel fiber material and decontamination tests with 
filters exposed to the WCF off-gas are being con- 
ducted, 


Laboratory and Pilot-Plant Investigations. A num- 
ber of laboratory and pilot-plant investigations are 
being carried out: WCF product storage, zirconium 
waste calcination, effect of NaK at high temperatures 
on stainless steels, and the testing of a valveless 
feed-control system. 

1. Effect of Storage on WCF Solids Product. The 
calcined solids from the WCF are stored underground 
in large stainless-steel bins.'! The center-line tem- 
perature of the stored calcine may be as high as 
700°C during storage due to the radioactive decay of 
the stored fission products. Constituents of the waste 
that are volatile at this center-line temperature mi- 
grate to a cooler zone where they condense. This re- 
sults in a change in the distribution of the volatile 
constituents during the period of high heat generation. 

Long-term (up to 1 year) heating experiments’’"’!° 
were conducted in which heated 2-g samples of 
actual WCF calcine were maintained at three tem- 
peratures: 525, 700, and 800°C. The results of these 
experiments were as follows: 

(a) Essentially no cesium volatilizes at 700°C or 
lower, but significant amounts volatilize at 800°C or 
higher. 

(b) Appreciable quantities of ruthenium volatilize 
at temperatures as low as 525°C. 

(c) The originally amorphous alumina converts to 
either gamma or alpha crystalline alumina, depend- 
ing upon the temperature. Gamma alumina is the 
only crystalline form at 700°C or lower, whereas 
alpha alumina is the predominant crystalline form 
at 800°C or higher. 

(d) No strontium or cerium is volatilized during 
storage at 800°C or lower. 


All the results obtained thus far indicate that no 
radionuclides will escape from the solids stored at 
high temperatures; instead, movement of the heat 
sources due to radionuclide migration may result 
in a stable situation in which the maximum tempera- 
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ture of the solids is lower thanthat without migration. 
Results of future studies planned to determine the 
interrelation between radionuclide migration and the 
resulting temperature distribution may show the 
feasibility of using even larger, lower-cost bins for 
storing calcined solids than the bins in the second 
storage facility'* currently being filled atICPP, 

2. Calcination of Zirconium Wastes. Wastes con- 
taining zirconium fluoride will be generated as a 
result of processing reactor fuels containing zirco- 
nium as one of the constituents. Studies*-!! are in 
progress to determine the possibility of calcining 
zirconium fluoride wastes near the end of the cur- 
rent WCF processing campaign. Tests are being 
conducted in the 12-in.-diameter pilot-plant calciner 
with a simulated zirconium fluoride waste solution 
being used as feed. Tests are also being carried out 
in a corrosion—erosion test loop to determine the 
effects of solids from the calcination of the zirco- 
nium fluoride wastes in the off-gas scrub solution. 

3. Performance of Stainless Steels in NaK at High 
Temperatures. A service test evaluation was con- 
ducted on a welded assembly* containing 304, 316, 
and 348 stainless steel.'® The assembly was inserted 
as a section of the NaK transmission returnline from 
the pilot-plant calciner and was exposed to NaK for 
14,200 hr without failure at temperatures from 
slightly below 1000 to 1400°F. There was no gen- 
eral corrosion damage on the materials. Appreci- 
able carbide precipitation at the grain boundaries 
and sigma phase (a hard intermetallic compound 
of the general formula Fe,,Cr,) were observed in the 
304 and 348 stainless steel; only carbide precipitation 
was noted in the 316 stainless steel. The appearance 
of the sigma phase, which is probably the most signif- 
icant deleterious effect, does not appear to affect the 
serviceability of the steel so long as the sigma-phase 
particles and the concentration are small.j These 
observations, together with evaluations of physical- 
property tests, indicate that these three stainless 
steels appear to be equally suitable from a corrosion 
standpoint for service with NaK in the temperature 
range and for the exposure time reported. 

4. Valveless Feed-Control System. Because of 
the plugging and erosion of valves in the WCF feed 
system, a valveless, nonrestrictive control system 
is being sought for future application£—10,12. Two ap- 
proaches were tested: air injection and a vortex- 
flow restrictor. Air injection alone is not satisfactory 
because it produces slugging flow. Results of small- 
scale tests of an air-operated vortex-flow restrictor 





*The assembly (nominal diameter of 2 in.) consisted 
of a schedule-40 348 stainless-steel elbow, one end of 
which was welded to a schedule-40 316 stainless-steel 
nozzle and the other end to a schedule-80 304 stainless- 
steel pipe. 

tIn applications where sigma-phase growth is undesir- 
able, the author suggests the use of a fully austenitic 
stainless steel with low carbon content. 


have been very encouraging. Additional tests will be 
made in a full-scale mock-up of the WCF feed system. 


Selected Current Literature 
FEED-INJECTION STUDIES 


An important feature in the operation of the WCF 
calciner is the feed injection of waste solutions into 
the heated, fluidized bed with pneumatic atomizing 
nozzles. In a recent article, Legler,'® Idaho Nuclear 
Corporation, reports on the feed-injection studies 
made in connection with the development of the WCF 
before the first waste-processing campaign. Of sev- 
eral schemes tested for injecting feed solution into 
the bed, best results were obtained with a two-fluid 
pneumatic atomizing nozzle unconventionally sub- 
merged in a heated, fluidized bed. Three such noz- 
zles (two-fluid, external-mixing atomizing nozzles 
of conventional design) were installed in the WCF 
calciner and successfully used in injecting plant 
wastes (first processing campaign) into the calciner. 
Their operation was considered very satisfactory. The 
author believes that this and other types of conven- 
tional pneumatic atomizing nozzles may be suitable 
for injecting feed into other fluid-bed processes where 
conversion of liquids to granular solids is desired. 


ASPHALT IMMOBILIZATION OF INTERMEDIATE-LEVEL- 
ACTIVITY WASTES 

Oak Ridge National Laboratory (ORNL) is devel- 
oping a process for the asphalt-immobilization of 
intermediate-level-activity wastes (ILWw). 20 Two re- 
ports2!»22 related to the development of this process 
have become available. The first report”! summarizes 
the laboratory studies conducted at ORNL on this 
process. It covers laboratory work concerned with 
the determination of process stoichiometry and pa- 
rameters and with the testing of asphalt? products 
for leaching rate of various elements, radiation 
stability, and safety evaluation. The other report?” 
presents the results of safety tests carried out by 
a commercial laboratory, '® under subcontract to 
ORNL, on asphalt products. Samples of pure asphalt 
(plus emulsifying agent) and asphalt products con- 
taining 20, 40, and 60 wt.% salts$ were subjected 
to various laboratory tests, including those of sta- 
bility (drop weight, autoignition, and shock sensi- 
tivity) and flame propagation. Included as part of 
the testing program were large-scale burning tests 
of 30-gal drums of asphalt products. In the laboratory 





tThe asphalts used in these studies contain nominally 
63 wt.% asphalt, 35 wt.% water, and 2 wt.% emulsifying 
agent. 

§In the process proposed by ORNL, the ILW is first 
concentrated in a waste evaporator. This concentrate?! is 
principally a nitrate solution of sodium, potassium, alumi- 
num, and iron and often contains sulfate, chloride, fluoride, 
and phosphate as well. 
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tests all materials tested proved to be insensitive to 
both impact and shock. However, it was not possible 
to establish the autoignition temperature and burning 
rate for any of the materials. In the large-scale 
burning tests, the salt—asphalt mixtures burned 
vigorously; they did not explode or show a significant 
increase in burning rate over that for the pure 
asphalt. 

The Nuclear Safety Information Center at ORNL 
has issued an indexed bibliography”® containing ab- 
stracts of articles dealing with the following subjects: 
shipping containers, shipping regulations, criticality 
safety as related to shipping and handling, transporta- 
tion accidents, and all other items dealing with safety 
during the transportation and/or handling of radioac- 
tive materials. 
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Neither the United States, nor the Commission, nor any person acting on behalf of the Com- 
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use of any information, apparatus, method, or process disclosed in this journal may not in- 
fringe privately owned rights; or 

B. Assumes any liabilities with respect to the use of, or for damages resulting from the 
use of any information, apparatus, method, or process disclosed in this journal. 


As used in the above, ‘‘person acting on behalf of the Commission’’ includes any employee 
or contractor of the Commission, or employee of such contractor, to the extent that such em- 
ployee or contractor of the Commission, or employee of such contractor prepares, dissemi- 
nates, or provides access to, any information pursuant to his employment or contract with the 
Commission, or his employment with such contractor. 
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